(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



,19,W„HdI.«™^P2^0r8..-«.» ■■lllllll Ilinillllllllllllll 

(43) International Publication Date (10) International Publication Number 

17 July 2003 (17.07.2003) pCT WO 03/057025 A2 



(51) laternational Patent Classification^: 

5/103, 5/04, 5/0452, G06F 17/00 



(74) Agents: LUZZATTO, Kfir ct al.; P.O. Box 5352, 84152 
Beer-Sheva (IL). 



(21) International Application Number: PCT/IL03/00005 

(22) International Filing Date: 2 January 2003 (02.01.2003) 



(25) Filing Language: 

(26) Publication Language: 



(81) Designated States (national): AH, AG, AL, AM, AT, AU, 

AZ, BA, BB, BG, BR, BY, BZ, CA, CH, CN, CO, CR, CU, 
CZ, DC, DK, DM, DZ, IiC, ED, ES, VI, GB, GD, GE, GH, 
GM, HR, HU, ID. IL, IN, IS, JP, Kti, KG, KP, KR, KZ, LC, 
LK, LR, LS, LT, LU, LY MA, MD, MG, MK, MN, MW, 
MX, MZ, NO, NZ, CM, PH, PL, PT, RO, RU, SO, SD, SE, 
SG, SK, SL, TJ, TM, TN, TR, IT, TZ, UA, UG, US, UZ, 
VC, VN, YU, ZA, ZM, ZW. 



147502 7 January 2002 (07.01 .2002) IL 

(71) Applicant (for all designated States except US): 
WIDEMED LTD, [HTIL]; P.O. Box 3002, Omer In- 
dustrial Park, 84965 Omer (IL). 

I (72) Inventors; and 

I (75) Inventors/AppUcants (for US only): GEVA, Amir 

i [HJTLi; 134 Hashalom Street, RO. Box 1589, 85025 

j Meitar (n.). TODROS, Kobi [nvnj; 79/2 Atad Street, 

i 84847 Beer Sheva (IL). PRESSMAN, Assaf [IL/ILl; 8 

■ Pnimt Street, 85538 Lehavim (IL). LEVI, Baruch [IL/IL]; 

; 2.3/28 Hativa Shmone Street, 842.'i2 Beer Sheva (H,). 

j MERGl, Aryeh [IL/IL]; 10 Rakefet Street, 45206 Hod 
(IL). 



(84) Designated States (regional): ARIPO patent (GH, GM, 
KR, LS, MW, MZ, SD, SL, SZ, TZ, UG, ZM, ZW), 
Eurasian patent (AM, AZ, BY, KG, KZ, MD, RU, TJ, TM), 
European patent (AT, BE, BG, CH, CY, CZ, DE, DK, EE, 
ES, FI, FR, GB, GR, IE, IT, LU, MC, NI., PT, SE, SI, SK, 
TR), OAPI patent (BP, BJ, CP, CG, CI, CM, GA, GN, GQ, 
GW, ML, MR, NE, SN, TD, TG). 

PubUshed: 

— without international search report and to be republished 

upon receipt of that report 

For nvo-leiwr codes and other abbreviations, refer to the "Guid- 
ance Notes on Codes and A bbreviations " appearing at the begin- 
ning of each regular issue of the PCT Gazette. 



(54) Title: SELF-ADAFITVE SYSTEM, FOR THE ANALYSIS OF BIOMEDICAL SIGNALS OF A PATIENT 

(57) Abstract: Method and system for monitoring normal/abnormal physiological events in patients by analyzing their biomedical 
signals, such as EEG signals, EOG signals, EMG signals, ECG signals and Respiratory signals. Corresponding raw data streams are 
aggregated, where each of which representing a different type of biomedical signals of the patient. Raw data streams are adaplively 
segmented and characterized for each type of biomedical signal. An attribute domain is determined, in which each segment is rcp- 
resented by a point corresponding to the attribute values of the segment. A set of clusters is generated in the attribute domain for 

^ each type of biomedical signal or for any combination thereof, each cluster consisting of a combination of points determined by their 
relative location to other points by assigning a set of property values to each point such that each property value corresponds to the 

2 degree of association of the point with one of the clusters. The combination of points included in each cluster is varied according to 
their most updated jffoperly value and by including points derived from adaptive segmentation of further aggregated raw data streams 

Q of the patient. At least one updated cluster is associated with a normal/abnormal physiological state of the patient by using former 

^ knowledge, regarding normal/abnormal physiological states of the patient and/or of a reference group of patients, that is represented 

^ as reference clusters in the domain. 
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SELF-ADAPTIVE SYSTEM. FOR THE ANALYSIS 
OF BIOMEDICAL SIGNALS OF A PATIENT 



Field of the Invention 

The present invention relates to the field of biomedical analyzing systems. More 
particularly, the present invention relates to a method and system for allowing 
sampling, processing and analyzing a plurality of aggregated (noisy) biomedical 
signals, such as Electrocardiogram (ECG), ElectroEncephalogram (EEG), 
Electro-Oculogram (EOG), Electromyogram (EMG), Oximetry and Respiratory 
signals, and for introducing to a clinician essentially noiseless and artifact-free 
biomedical event-wise and time-wise selections of analyzed data related to 
specific patient being monitored. The present invention also relates to a system 
that allows integrated analysis of ECG, EEG, EOG and respiratory signals 
while said system automatically adapting itself to the inherent biomedical 
parameters of said monitored patient. The present invention also relates to a 
Remote Medical Care system that allows remote monitoring of a person, by 
employing "Smart Compression" process on relevant physiological data and 
transferring the compressed data to remote end, for completing the evaluation 
process of the compressed data. 

Backsaround of the Invention 

Biomedical signals "suffer" from common problem of poor signal-to-noise ratio 
(SNR). This problem becomes acute whenever it comes to measuring 
ElecroEncyphalogram (EEG) and ElecroMyogram (EMG) signals, since the 
magnitude of such signals is generally in the fiVolt range. The signals picked- 
up by corresponding detectors, which are attached to different portions of the 
patient body, are forwarded to the analyzing system by wires, which may add 
significant noise to the already weak signal. In addition, the characteristics of 
the wires may change according to the ambient conditions, a problem that must 
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be overcome by utilizing special hardware and/or software resources. Another 
problem arises whenever the patient changes his relative position. Under such 
circumstances, the analyzing system may interpret false signals as valid 
pathological events. Therefore, it wotdd be advantageous to have an analyzing 
system that is capable of recognizing and disregarding such false events. 

In some conventional systems, wherein several physiological parameters of a 
patient are examined (e.g. EEG, ECG and EMG and respiratory), each 
parameter is analyzed regardless of other physiological parameters. However, 
in some cases, several parameters could be strongly interrelated. For example, 
an 'irregular' heartbeat rate (HBR) of a patient might be interpreted by a 
chnician as pathological HBR if analyzed regardless of other factors. However, 
it might be interpreted as proper HBR if correlated with brain waves indicating 
to the analyzing system that the patient is dreaming. Currently, in such 
systems a clinician has to 'manually' correlate the various raw measurements, 
which are related to different physiological aspects, and reach conclusions 
therejfiTom. 

Therefore, it would be advantageous to have an analysis system, which would 
allow measuring several types of physiological signals, and coiild introduce to 
the clinician already-correlated data, i.e., a data that is already correlated with 
(other) relevant physiological parameters so as to eliminate artifacts and other 
types of noises/interference. Such a system would introduce to the chnician 
relevant, more accurate and well-characterized selected physiological events, 
and save the clinician analyzing/interpretation time. For example, it would be 
advantageous to combine abdominal effort. Oxygen saturation, patient's, 
movement or arousals, and sleep stages ia order to reach a decision regarding 
apneas. 

Different patients are hkely to introduce different 'normal' physiological 
activities due to the variance of, e.g., heartbeat morphology (see Fig. la and lb). 
Therefore, it would be advantageous to have a system that would be capable of 
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automatically adapting itself to the specific heartbeat morphology of a person 
being monitored, because such system would allow more acc\irate detection of 
physiological abnormalities. 

Other conventional systems, generally called 'semi-automatic' systems, are 
capable of sampHng, processing and analyzing electro-physiological signals. 
However, these systems configuration allows them to only analyze physiological 
signals between two 'rigid' boundaries, in time and/or in magnitude, which 
must be predetermined by the clinician. The two predetermined boundaries are 
selected on a random base, and, therefore, processing and analyzing the raw 
data contained within these boundaries usually does not yield satisfactory 
results. In some cases, the results could be even meaningless. Therefore, the 
clinician must iterate a sequence of operations, on a 'trial and error' basis, until 
meaningful results are obtained. Therefore, it would be advantageous to have a 
system that would be capable of 'self-tuning' to the required data portion, 
without needing to depend on specific or predetermined boundaries of any type, 
between which the raw data is to be processed and analyzed. 

There are several common diseases associated with heart conditions, which are 
characterized by the corresponding irregiilarities in the heart behavior between 
distinguishable points in a heartbeat cycle (commonly known as points 'P', 'Q', 
'R', 'S' and 'T'). Figs. 1(a) and 1(b) illustrate typical heartbeats, in which the 
'P', 'Q', 'R', 'S' and 'T' points are marked. Conventional systems utilize 
processes of rather limited flexibility (e.g. using 'rigid' predetermined 
boundaries). Therefore, they offer solutions that focus on limited aspects of 
heart analysis. For example, one system focuses on analyzing irregularities 
between points 'S' and 'T', another system on irregularities between points 'P' 
and 'R' etc. It would be advantageous to have a system that would be capable of 
isolating and classifying essentially every type of irregidarity associated with a 
heart condition, regardless of the deformation (see Figs. 1(b), 1(d) and 1(e)) and 
relative location of such irregularity. 



wo 03/057025 PCT/n.03/00005 
-4- 

Conventional systems usually display measurement results graphically on 
paper stripe or on electronic display, after which they are reviewed and 
estimated by a chnician or therapist. In some cases, these results are stored in 
an electronic storage media. Sometimes, the physiological activity of a patient is 
monitored for several hours but the clinician/therapist is interested only in 
small portions of the recorded data, which are associated with abnormal 
physiological activity of the patient. In conventional systems, the rest, and 
sometimes most, of the recorded data is useless because it is either too 
corrupted or it does not contribute anything to the assessment of the patient 
physiological condition, and is, therefore, dumped away (whenever a paper 
stripe is used), or a data storage place is unnecessarily consumed (whenever an 
electronic storage media is used), which leads to wasted resources. Other 
applications require collecting physiological signals of specific patient for 
several hours in order to allow good assessment of the patient physiological 
condition. In this case, the clinician must review and analyze an enormous 
amount of data, which is usually a burdensome process. The analysis lasts a 
long period of time and changes with technician tiredness. It would be 
advantageous to have a system, which is capable of performing the analysis 
with the same reUabihty. 

In connection with sleep disorders, most current sleep (polysomnographic) 
studies are carried out in Sleep Centers, where patients need to spend at least 
one night in a hospital or private laboratory. The fact that patients need to 
sleep in the laboratory facilities presents significant difficulties. First, with the 
growing awareness to sleep disorders, waiting lists are growing, and delays of 
months and years are common. Second, the sleep study is expensive and thus 
the availability of sleep monitoring is more difficult. But more important fact is 
that the patient's sleep is affected whenever studied in a non-familiar 
laboratory setting. For example, sleep efficiency and patterns are badly affected 
due to the phenomena generally known as the "first night" effect. Therefore, in 
order to achieve a genuine sleep behavior of a patient, it is necessary to study 
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the patient (especially for clinical research) for at least two consecutive nights, a 
fact that increases the expenses. 



The advantages of conducting sleep studies at the patient home include 
increased flexibility of scheduling, less disruption to the routine of the child and 
family, and the ability to study the patient, especially a child, in his, or her, 
natural environment. While this has not yet been documented, it seems that 
home studies will be more cost-efficient than laboratory polysomnography, as a 
hospital admission is not required. Moreover, because such studies usually do 
not require the attendance of technicians, technician time is substantially 
decreased, or it could even be evaded. These factors become even more 
important when more than one night of recording is reqiiired. 

More and more Sleep Centers are trying to screen their patients by performing 
partial ambulatory sleep studies at home. Although this is a growing field in 
sleep medicine, lacking automatic signal processing and monitoring greatly 
complicates the amattended tests. 

A home video and partial cardiovascular electrode are successfully employed in 
Montreal, using a customized computer program. The possibUity of using a 
simple, non-invasive data acquisition at home, monitored via -a 
computer/Internet, on line, may facihtate home sleep monitoring for millions of 
patients. A patient will be able to remain in his own bedroom, for as long as 
required, while being continuously monitored and the data inspected on-line by 
a technician, enabUng him to handle a problematic signal from the center or 
even call the patient at his home, to prevent waste of time due to disrupted data 
acquisition. In addition, the data will be analyzed on-hne, in real time and 
automatically by the computer, thereby saving a lot of technician working hours 
and making the results readily available for both patient and his clinician. 

Conventional systems produced by companies hke Compumedics, Flaga, 
Respironoics, etc., are capable of automatically analyzing the mentioned 
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signals. However, in order to achieve reliable results, an array of thresholds 
should be provided and adjusted manually in order for these systems to 
overcome the problem of different patients having different inherent signal 
variations and, thereby, to manually converge to the patient self-parameters. 

There are major three frequency bands which characterize the EEG signal 
during sleep. The statistical characteristics of the frequency content of each 
group statistically differ from one person to another. Therefore, it is essential 
that a sleep monitoring system would have an adaptive mechanism, which 
would allow it to reach correct decisions regarding the interpretation of each 
frequency group on individual basis. 

All of the methods described above have not yet provided satisfactory solutions 
to the problem of automatically adapting and optimizing the biomedical 
analyzing system to a patient being monitored. Additionally, the prior art 
methods have not yet provided satisfactory solutions .to the problem of 
integrating analysis of ECG, EMG, EEG and respiratory signals. 

It is an object of the present invention to provide a system that is capable of 
automatically adapting/optimizing itseffto the patient being monitored. 

It is another object of the present invention to provide a method and system for 
allowing integrating analysis of ECG, EEG and respiratory signals. 

It is still another object of the present invention to provide a method and system 
for allowing extracting from same electro-physiological signal several data 
types, each of which being related to a different type of physiological activity. 

It is a further object of the present invention to provide a method for allowing 
introducing to a clinician, according to his request, only selected biomedical 
data. 
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It is yet another object of the present invention to provide a method and system 
that allows efficient transmission of selected biological data to remote locations. 

It is still another object of the present invention to provide a system that is 
capable of identifying sleep stages of a patient. 

It is still a further object of the present invention to provide a method for 
allowing introducing to a clinician, upon request, data representing the 
physiological state of a patient, derived from his biomedical data. 

Other objects and advantages of the invention will become apparent as the 
description proceeds. 

Summary of the invention 

The present invention is directed to a method for monitoring normal/abnormal 
physiological events in patients by analyzing his biomedical signals, such as 
EEG signals, EOG signals, EMG signals, EGG signals and Respiratory signals. 
For each patient, one or more raw data streams are aggregated, where each of 
which representing a different type of biomedical signals of the patient, in an 
accessible data storage. Segmentation of the raw data streams is performed 
adaptively to the patient for each type of biomedical signal, and individual 
attributes that are represented by attribute values are assigned to each 
segment. An attribute domain, in which each segment is represented by a point 
that corresponds to the attribute values of the segment, is determined. A set of 
clusters is generated in the attribute domain for each type of biomedical signal 
or for any combination thereof, each cluster consisting of a combination of 
points determined by their relative location to other points by assigning a set of 
property values to each point such that each property value corresponds to the 
degree of association of the point with one of the clusters. The probability of the 
order of appearance of each point, in time, is obtained according to its property 
value and the property value of each point in each cluster is updated using the 
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obtained probability. Each cluster is updated by repeating this process while in 
each time, the combination of points included in each cluster is varied according 
to their most updated property value and by including points derived from 
adaptive segmentation of further aggregated raw data streams of the patient. 
At least one updated cluster is associated with a normal/abnormal physiological 
state of the patient by using former knowledge, regarding normal/abnormal 
physiological states of the patient and/or of a reference group of patients, that is 
represented as reference clusters in the domain. The patient is individually 
characterized by identifying normal/abnormal physiological states, associated 
with one or more updated clusters, and obtaining the probability of a change 
between normal/abnormal physiological states using the probability of the order 
of appearance. The order of appearance of the normal/abnormal physiological 
states are utilized for predicting normal, and/or abnormal, physiological states 
of the specific patient. 

EEG signals are segmented and the segments are clustered (the segmentation 
and clustering processes are described herein below) and the clusters are 
associated with, e.g., three different frequency groups, each of which being 
related to different corresponding sleep stages. Each cluster of EEG signal 
segments is associated with at least one corresponding frequency group. 
Clustering of segments of the EEG signal may be carried out according to a 
different number of (i.e., different than three) frequency groups. 

Relative activity of each frequency group is calculated within a moving time- 
window, or time frame, of 30 seconds. The time-window is adaptively changed 
(i.e., in accordance with temporal and/or frequency attributes of the EEG 
signal). 

The frequency groups may be for example: high-frequency group, which 
represents Wake and REM stages, mixed-frequency group, which represents 
sleep stages 1 and 2, and low-frequency group, which represents sleep stages 3 
and 4. Accordingly, each time window 'contributes' three average values that 
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represent, in that time window, the relative activity of each frequency group, 
with respect to the activity of the other frequency groups. A graph, or the 
behavior, of each one of the frequency groups may be formed, or obtained, by 
utilizing corresponding consequent average v£dues of the respective frequency 
group. 

The present invention offers a novel method for displasdng the aforesaid 
relative activities, or graphs, i.e., the relative activities of the high frequency 
group, mixed frequency group and low frequency group, are displayed as a 
function of time, in the form of two-dimensional color-coded graph. In the graph, 
each one of the frequency groups is displayed by a different color, for allowing 
an observer (i.e., a phs'^sician) distinguishing between the different frequency 
groups. The color-coded graph provides to an observer detailed information 
related to the dynamics of a sleep of a monitored person, as well the 'purity' of 
each sleep stage, which are essential for diagnosing sleep disorders. 

Indications related to abnormal physiological events in the patient are obtained 
by displaying changes and/or different stages in normal physiological states of 
the patient as a function of time. At least a portion of the aggregated raw data 
streams may be stored in accessible data storage. Data segment that 
corresponds to an identified disorder may be retrieved from the accessible data 
storage. 

Preferably, individual data being a group of updated clusters and/or indications 
of normal/abnormal physiological events in each patient represented by the 
group are generated. Each individual data is divided into packets that are 
transmitted to a remote site over a data network. The data packets are received 

in the remote site, processed by a computerized system and a representation of 
normal/abnormal physiological events is displayed as a function of time and/or 
of changes in normal physiological states and/or different stages of 
normal/abnormal physiological states of each patient at the remote site. 
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A request for one or more segments of the raw data streams that correspond to 
abnormal physiological events may be transmitted from the remote site to the 
data storage over the data network. The requested segments are retrieved from 
the data storage and displayed at the remote site. Instructions for a medical 
treatment to be given to the patient for terminating the abnormal physiological 
events may be transmitted over the data network from the remote site to a local 
computerized system located at the vicinity of the patient. Data packets that 
correspond to an individual group of data clusters that corresponds to 
physiological events after applying the medical treatment to the patient are 
received at the remote site, and processed by the computerized system. A 
representation of physiological events, obtained in response to the medical 
treatment, is displayed as a function of time and/or of changes in normal 
physiological states and/or different stages of normal physiological states of 
each patient. This process is repeated until abnormal physiological events are 
eliminated or substantially reduced. 

Raw data streams may be received from one or more dry electrodes attached to 
the skin of the patient and movement detectors may be utilized for identifying 
time periods during which segments, or portions, of raw data streams are 
associated with (i.e., affected by) patient movements. 

The abnormal physiological events may comprise sleep apnea, cardiac 
arrhythmia, Ventricular Tachycardia, Ventricular Fibrillation/Flutter and 
epilepsy. 

In one aspect, the raw data streams may be preprocessed by inputting the raw 
data streams to a cellular device having processing means that are capable of 
identifying irregular physiological events of the patient. A characteristic trace 
that corresponds to regular physiological events of the patient is generated by 
the cellular device that identifies selected raw data segments that comprise 
irregular physiological events and transmits the selected raw data segments 
over a cellular network to another computational device. The selected raw data 
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segments are further processed and the processed raw data is represented by 
data clusters, each of which consisting of a set of xmique properties that 
correspond to a different type of biomedical signal and to a different normal 
physiological state of the patient, by another computational device. 

Preferably, the abnormal physiological events and/or the individual attributes 
of each patient are represented by an indiAddual group of updated data clusters, 
which are divided into data packets that are transmitted to a remote site over a 
data network. The data packets are received in the remote site, processed by a 
computerized system and a representation of abnormal physiological events is 
displayed as a function of time and/or of changes in normal physiological states 
and /or different stages of normal physiological states of each patient at the 
remote site. 

Preferably, irregular physiological events are identified for each patient, by 
another computational device. Segments of raw data that corresponds to the 
irregular physiological events are retrieved and a data set that comprises the 
retrieved segments and parameters indicating the irregularity of the irregular 
physiological events within the retrieved segments is generated. The data set is 
transmitted to a cellular device having processing means which are capa,ble of 
displaying the parameters and the retrieved segments, over a cellular data 
network. The parameters and/or the retrieved segments are then displayed on 
the display of the cellular device. 

In one aspect, alert indications representing abnormal physiological events are 
generated in real-time and automatically transmitted to a physician over a data 
network, for determining/providing the required medical treatment. 

The operation of medical apparatus (such as, controllable drug dosage devices, 
controllable gas/air delivery devices, and CPAPs) used for providing medical 
treatment to a patient being in an abnormal state may be dynamically 
controlled in response to identified abnormal physiological events in the patient. 
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Instructions for dynamicEilly controlling medical equipment, providing medical 
treatment to the patient, for terminating the abnormal physiological events 
may be transmitted from the remote site to a local computerized system located 
at the vicinity of the medical equipment, over the data network. 

The former knowledge may be capable of being continuously updated by 
manually modifying the characterization results, e.g., by a physician examining 
these resxdts. 

The present invention is also directed to a system for monitoring 
normal/abnormal physiological events in one or more patients by the analysis of 
biomedical signals of each patient, that comprises: 

a) One or more sensors, attached or allocated to the patient, being capable 
of generating and transmitting biomedical signals; 

b) a data storage for aggregating transmitted raw data streams, each of 
which representing a different type of biomedical signals of the patient; 
and 

c) processing means, being in data communication with the sensors and/or 
the data storage, for performing adaptive segmentation of the raw data 
streams of the patient and assigning, to each segment, individual 
attributes being represented by attribute values determining an attribute 
domain, in which each segment being represented by a point, for 
generating a set of clusters in said attribute domain for each type of 
biomedical signal or for any combination thereof, such that each cluster 
consists of a combination of points determined by their relative location 
to other points, by assigning a set of property values to each point, and 
such that each property value corresponds to the degree of association of 
said point with one of the clusters, for obtaining the probability of the 
order of appearance of each point, in time, according to its property value 
and updating the property value of each point in each cluster using said 
probabihty, for updating each cluster by iteratively varjdng the 
combination of points included in each cluster according to their most 
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updated property value and by including points derived from adaptive 
segmentation of further aggregated raw data streams of said patient, for 
associating at least one updated cluster with a normal/abnormal 
physiological state of said patient by using former knowledge, regarding 
normal/abnormal physiological states of said patient and/or of a reference 
group of patients, that is represented as reference clusters in said domain 
and for individually characterizing said patient by identifying 
normal/abnormal physiological states, associated with one or more 
updated clusters, and obtaining the probability of a change between 
normal/abnormal physiological states using said probability of the order 
of appearance. 

The system may further comprise display means for displaying different stages 
of normal physiological states of the patient as a function of time. 

Preferably, the system further comprises: 

a) Means for dividing data clusters and/or the representation of their 
corresponding physiological states into data packets; 

b) means for transmitting the data packets of one or more patients to a 
remote site over a data network; 

c) means for receiving the data packets in the remote site; and 

d) means for processing the received data packets and displaying a 
representation of abnormal physiological events as a function of time 
and/or of changes in normal physiological states and /or different stages of 
normal physiological states of each patient at the remote site. 

The system may further comprise: 

a) means for transmitting a request for one or more segments of the raw 
data streams that correspond to abnormal physiological events from the 
remote site to the data storage over the data network; and 

b) means for retrieving the requested segments from the data storage and 
displaying the requested segments at the remote site. 
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One or more dry electrodes having movement detectors for generating and 
forwarding raw data streams may be attached to the skin of the patient. 

The system may further comprise: 

a) a computational device having processing means for identifying abnormal 
physiological events in a patient, and being capable of generating a data 
set that comprises retrieved segments and parameters indicating the 
irregularity of the irregular physiological events within the retrieved 
segments and of transmitting the data set; and 

b) a cellular device, having processing means being capable of receiving the 
data set over a cellular data network, and a display for displaying the 
parameters and the retrieved segments. 

The following terms will be hereinafter used: 

1. 'Gateway' - a Gateway may be a personal computer (PC), to which an IP 
address is pre-assigned, having a corresponding software package, for 
allowing it to perform the following tasks: (a) reading, or accepting, raw data 
from corresponding biomedical monitoring devices capable of recording 

biomedical signals of patients connected thereto; (b) forwarding to patient 
treatment devices digital, or analog, control signals, for affecting the 
physiological status of patients connected thereto (by, e.g., causing the 
treatment device to administer drugs into the patient; (c) formatting the raw 
data iato corresponding format, which could be a proprietary format or any 
known suitable format, for allowing communicating the formatted data to a 
(remote) Lab Server. 

2. 'User Workstation' - a user workstation is a computerized system (e.g., a 
personal computer), to which an IP address is assigned, allowing users, or 
clients, such as physicians and other medical and clinical staff, forwarding 
requests to a Lab server, or, optionally, to a central server, via the Internet 



wo 03/057025 PCT/n.03/00005 
-15- 

HTTP protocol, for retrieving and displaying preferred data related to 
preferred aspects of physiological condition of patients. 
3. 'Lab Server' - the Lab Server is a computerized apparatus which includes at 
least a data base, data processing unit and URL (Uniform Resource 
Locator) address for allowing plurality of distributed gateways and user 
workstations an access via the Internet. The Lab Server tasks are at least: 
(a) identifying users by utilizing a table of users, and allowing users to 
view/change various portions of content, of medical, administrative and 
control files contained therein, according to the users' profiles (i.e., 
authorization level); (b) Receiving data/information firom Gateways, and 
classifying and storing the received data; (c) processing and analyzing the 
data/information according to the processes disclosed in the current 
invention, and storing the analyzed data and related statistics data; (d) 
Generating, from data resulting firom the analyzed data, HTML/XML files, 
according to user's request; and (e) forwarding specific data (optionally 
encrypted) in the form of packets, HTTP, HTML or XML, according to the 
user's request (i.e., the request may refer to study name, time interval, 
preferred aspects or events of physiological condition of a patient, biomedical 
channel (e.g., heart beat channel, etc.); (f) forwarding to gateways requests 
for specific data, and/or control signals to the monitoring/treatment devices. 
The requests and control signals may be initiated by a user, by utilizing the 
user's workstation. The control signals to a monitoring device could relate 
to, e.g., the initialization parameters of the monitoring device. Preferably, a 
Lab Server contains a database, which includes data related to at least: 
identification details (IDs) of all the monitoring and treatment devices 
connected thereto; clinical records per patient; patients' personal details 
(e.g., name, ID, address, etc.); physicians of patients; statistics of clinical 
disorders and irregularities (i.e., 'events'); and authentication details, for 
allowing an access to the system only to authorized personnel. Of course, 
other relevant data could be stored in the database, and the database could 
be updated whenever required. 
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4. 'Central server' is a computerized apparatus which, includes at least a 
central database, URL and application software, the tasks of which are 
described hereinbelow. Each one of the Lab Servers and, optionally, user 
workstations, may communicate with the (remote) central server, the tasks 
of which are: (a) Defining authorization profiles of user workstations and 
users utilizing the user workstations, and the relation between the users. 
Lab Servers and patients; (b) identifying users and the users' relation to 
specific Lab Server, by utilizing users', Lab Servers and patients' clinical 
and personal data reference table(s); and (c) routing or relaying users to 
specific Lab Server, from where requested data is retrieved and forwarded to 
the respective (requesting) users. 

According to a first most preferred embodiment of the present invention, there 
is provided a 'one-Lab-Server' system for allowing a plurality of users, each of 
which coxdd be connected to the system firom a different location, retrieving and 
displasdng, possibly simultaneously, data related to preferred aspects of 
physiological and medical condition of a plurality of patients, each of which 
covid be connected to the system from a different location, not necessarily the 
same locations as the location of the users, and affecting the physiological 
status, or condition, of patients connected to the system, comprising: 

a) at least one gatewas'', to which a plurality of monitoring, and, 
optionally, treatment devices, could be connected, each of which having a 
imique ID code. Each gateway is capable of preprocessing multiple raw 
data streams that represent the biomedical signals of patients connected to 
respective monitoring device, and controlling the operation of each one of 
the monitoring and treatment devices, by utilizing the aforesaid ID codes; 

b) at least one user workstation, for allowing a user (normally a 
physician) obtaining medical information of specific patient(s) and affecting 
the medical, or physiological, condition of the patient(s); and 

c) Lab Server. The main task of the Lab Server is processing and storing 
the data forwarded form the gateway(s) and forwarding requested data to 
corresponding User Workstation(s). Another task of the Lab server is 
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controlling the monitoring apparatuses that axe connected to the gateway, 
by sending control signals generated by the Lab server, and/or by a User 
workstation. The Lab Server is capable of communicating with the at least 
one gateway and with the at the least one user workstation, via a data 
network, such as the Internet. To the Lab Server is forwarded, via the data 
network, formatted raw, or preprocessed, data that represents the 
physiological condition of the at least one patient. The Lab Server is 
programmed to process and analyze the formatted raw, or preprocessed, 
data according to the novel method disclosed in the current invention, and 
the at least one user workstation is capable of forwarding to the Lab 
Server, via the data network, a request for retrieving, from the Lab Server, 
and displaying, on the user workstation, processed and analyzed data 
related to specific aspects of the physiological condition of a wanted 
patient. The specific aspects are normally defined by the users. Users may 
also request for data related to specific the statistics of various preferred 
physiological aspects. 

The present invention also provides a method for operating the one-Lab-Server 
system. Preferably, the method comprises: (a) login into a corresponding Lab 
Server using the user's name and password (encrypted communication is an 
option); (b) sending a request, from the user workstation to the lab server, for 
data related to specific physiological aspects of specific patient, and browsing, 
by a user (e.g., physician), through the pages of information that the user is 
authorized to view; (c) receiving the data and utilizing Graphic User Interface 
(GUI) for displaying recorded information (that was recorded, e.g., overnight) of 
bio-signals, such as hypnogram, heart rate, etc., every predefined time 
intervals, for example, displaying data each 30 seconds. The GUI allows also 
displaying zoom pictures of charts/graphs; (d) whenever required, defining 
changes that should be made (e.g., to compensate for known errors) to specific 
(physiological) data, and updating the Lab Server with these changes; and (e) 
writing textuEd information and comments, including medical interpretation, 
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and relating these textual information and comments to corresponding 
physiological data. 



According to another most preferred embodiment of the present invention, there 
is provided another type of system, (hereinafter referred to as a 'Mnlti- 
patientV'multi-physician' system), which comprises: 

a) at least two 'one-Lab-Server' systems, the configuration and functionahty of 
which Eire essentially similar to those of the one-Lab-Server system 
described above; and 

b) Central server. The main task of the central server is intermediating and 
coordinating between Lab Servers and user Workstations. The central 
server is capable of bi-directional communications with the Lab Server of 
each one of the at least two one-Lab-Server systems, via the data network, 
and the at least one user workstation of each one of the one-Lab-Server 
systems is capable of bi-directional communications with the respective Lab 
Server, or with the central server, or with both respective Lab Server and 
central server. The bi-directional communications includes forwarding, from 
any one of the user workstations to the central server, requests for data 
related to preferred aspects of the physiological condition of specific patient, 
regardless of the gateway to which the specific patient is or was connected, 
and, if the at least one request is forwarded directly to the central server, 
searching by the central server for the requested data in each one of the Lab 
Servers. Each one of the user workstations may forward a request for 
specific biomedical data to the Central server, to which the Central server 
responds by searching in the central database for a data that will be helpful 
in routing the user to the corresponding Lab Server (i.e., the Lab Server 
containing the requested biomedical data), from which the user may retrieve 
the required medical data. 



The central server functions to intermediate between all user workstations and 
all Lab Servers, and has URL and Internet network connection and has also a 
central database that contains data related to each one of the databases 
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residing within, or related to, every Lab Server, and in particiilar, the central 
database contains data related to users, their login information and 
authorization level information for each user; gateways parameters; and 
patients' personal and clinical data. The Central Server also functions to 
identify users, and to route or relay users to relevant Lab server(s), from which 
requested biomedical data is retrieved by the users. The Central server also 
functions as the general manager and coordiaator of the multi-patient/multi- 
physician system, thereby allowing automatic routing (connection) of users to 
corresponding Lab Server(s) by the Central Server. To each one of the Lab 
Servers coidd be connected, via a data network, a plurality of Gateways and 
user workstations, essentially iu the same manner as described in connection 
with the one-Lab-Server system. 

The identification of users by the Central Server may be performed by 
requesting and receiving corresponding information from an Authentication and 
Encryption Server. The Authentication server provides information for 
authentication of users and encryption keys for sectired transmission. 
Whenever a user with certain permission level logs-in (via the respective user 
workstation), the Central server forwards the 'login data' to the Authentication 
and Encryption server for authentication. The Authentication and Encryption 
server checks the permission level of the user, and allows the user to interact 
with the system (e.g., retrieving biomedical data and initiating control signals) 
according to this permission level. 

Of course, individual 'one-lab server' systems may be configured to be self- 
contained. 

Brief Description of the Drawings 



The above and other characteristics and advantages of the invention wUl be 
better understood through the following illustrative and non-limitative detailed 
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description of preferred embodiments thereof, with reference to the appended 
drawings, wherein: 

Figs, la to le (prior art) schematically illustrate typical normal 
heartbeats shape and two commonly known aberrant heartbeat shapes; 

Fig. 2 is a block diagram illustrating the system's adaptation to the 
Patient-Specific Universal Heartbeat (PSUHB), according to a preferred 
embodiment of the present invention; 

Fig. 3 is a block diagram illustrating the system's adaptation to the 
Patient-Specific Universal Pathological Heartbeat (PSUPHB), according to a 
preferred embodiment of the present invention; 

Fig. 4 schematically illustrates a more detailed system's adaptation to 
the Patient-Specific Universal Heartbeat (PSUHB), according to a preferred 
embodiment of the present invention; 

Fig. 5 schematically illustrates the heartbeat classification process, 
according to the preferred embodiment of the invention; 

Fig. 6 shows testing results of the first heartbeat classification process 
according to which viable heartbeats are classified, according to a preferred 
embodiment of the present invention; 

Fig. 7 schematically illustrates the pathological heartbeat first 
classification process, according to the preferred embodiment of the invention; 

Fig. 8 illustrates R-peak detection by employing the Wavelet process on 
ECG signal, according to a preferred embodiment of the present invention; 

Fig. 9 schematically illustrates the segmentation process of Normal 
heartbeats, according to a preferred embodiment of the present invention; 

Figs. 10a to lOe show test results that emphasize the various segments 
of Normal heartbeats, according to the preferred embodiment of the present 
invention; 

Figs. 11a to lib show test results that were obtained by monitoring four 
persons having different types of arrhj^hmias, according to the preferred 
embodiment of the present invention; 
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Fig. 12 schematically illustrates the process of detecting a large 
accidental peak that was superimposed on heartbeats, according to a preferred 
embodiment of the present invention; 

Fig. 13 schematically illustrates adaptive detection of low frequency 
movements that are superimposed on viable heartbeats, according to a 
preferred embodiment of the present invention; 

Fig. 14 schematically illustrates adaptive detection of high frequency 
EMG and noise signals that are superimposed on viable heartbeats, according 
to a preferred embodiment of the present invention; 

Fig. 15 illustrates a general block diagram of an automatic Sleep 
Staging System, according to a second embodiment of the invention; 

Fig. 16 schematically illustrates filtering of an EEG signal, according to 
a preferred embodiment of the present invention; 

Fig. 17 schematically illustrates adaptive segmentation process of an 
EEG signal, according to a preferred embodiment of the present invention; 

Fig. 18 schematically Ulustrates detection of frequency activities by 
utihzing a three-state HMM model (a) and a two-state HMM model (b), 
according to a preferred embodiment of the present invention; 

Fig. 19 (prior art) illustrates two quasi-stationary processes that were 
derived from an EEG signal; 

Fig. 20 is a general flow chart of describing the sequential segmentation 
process, according to a preferred embodiment of the present invention; 

Fig. 21a schematically iQustrates the principle of the initial boimdaries- 
detection of a quasi-stationary segment, a preferred embodiment of the present 
invention; 

Fig. 21b schematically illustrates the principle of optimization of a 
boundary location, according to a preferred embodiment of the present 
invention; 

Figs. 22a and 22b schematically Ulustrate envelope detection of 
Respiratory signal, according to the second embodiment of the invention; 

Fig. 23 illustrates utUization of parabola for signal trend evaluation, 
according to a preferred embodiment of the present invention; 
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Fig. 24a and 24b illustrate apnea events and corresponding decrease in 
Oxygen saturation level, according to a preferred embodiment of the present 
invention; 

Fig. 25 illustrates adjoining Down parabola and Up parabola to one 
parabola, according to a preferred embodiment of the present invention; 

Fig. 26 schematically illustrates exemplary combinations of different 
events, on which the Decision-Making Process (DMA) is employed, according to 
a preferred embodiment of the present invention; 

Fig. 27 (prior art) illustrates exemplary resulting apnea events 
according to a human scorer; 

Fig. 28 illustrates exemplary resulting apnea events, according to a 
preferred embodiment of the present invention; 

Fig. 29 schematically illustrates the general functionahty of the system, 
according to a preferred embodiment of the present invention; 

Fig. 30 schematically illustrates a Cellidar-phone-based Remote 
Monitoring System (CRMS), according to one embodiment of the present 
invention; 

Fig. 31 schematically illustrates utihzuig the Automatic Sleep Stage 
(ASS) system with combination with a Continuous Positive Airway Pressure 
(CPAP) system, according to a preferred embodiment of the present invention; 

Figs. 32a and 32b illustrate exemplary results of the adaptive 
segmentation of EEG and EOG signals, respectively, according to a preferred 
embodiment of the present invention; 

Fig. 33 illustrates an example for associating quasi-stationary segments 
with one of the three frequency groups, according to a preferred embodiment of 
the present invention; 

Fig. 34 illustrates detection of special events, according to a preferred 
embodiment of the present invention; 

Figs. 35a to 35d illustrate exemplary classification of special events, 
according to a preferred embodiment of the present invention; 

Fig. 36 illustrate exemplary REM events, according to a preferred 
embodiment of the present invention; 
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Fig. 37 illustrates exemplary clusters, according to a preferred 
embodiment of the present invention; 

Fig. 38 depicts the residting exemplary hypnogram, which describes the 
relevant sleep stages, according to a preferred embodiment of the present 
invention; 

Figs. 39a and 39b illustrate exemplary detection of low energy portions 
in the EMG signal, according to a preferred embodiment of the present 
invention; 

Figs. 40a to 40c illustrate typical Awakening, Arousal and Micro- 
Arousal stages, according to a preferred embodiment of the present invention; 

Fig. 41 is an exemplary 2-dimensional representation of relative activity 
of EEG frequency groups, according to a preferred embodiment of the present 
invention; 

Fig. 42 schematically illustrates basic 'one-Lab-Server' system, 
according to one most preferred embodiment of the present invention; and 

Fig. 43 schematically illustrates a system that comprises several 'one- 
Lab-Server' systems, according to another most preferred embodiment of the 
present invention. 

Detailed Description of Preferred Embodiments 

In one aspect, the present invention is directed to a biomedical- analyzing 
platform (hereinafter referred to as "TeleMedicine Platform" - TMP) that is fast 
responsive, patient-adaptive and having high quality ambxdatory remote 
monitoring of human biomedical signals, which are associated with cardiac 
status of a patient, his sleep patterns and respiratory behavior, by utUizing 
several known processing processes (i.e. HMM, Wavelet Transform process and 
Viterbi) and novel Adaptive Segmentation Processes (ASA) and fuzzy-logic 
based decision-making process. The present invention is also directed to a 
system that allows integrating analysis of a plurality of physiological signals. 
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General Introduction 

The present invention is characterized hy introducing a system that is patient- 
adaptive. According to the present invention, the system is capable of 
automatically adapting itself to the inherent cardiac, sleeping and respiratory 
status of individual patients, without having to manually match new sets of 
thresholds (on a 'trial and error' basis) for different patients being monitored. 
Therefore, the system optimizes itself to the persons being monitored on 
individual basis. 

The present invention is further characterized by the way it utilizes known 
processes. In connection with the ECG signal, the present invention utilizes the 
(known) Wavelet Transform Process (WTA) for adaptive segmentation of 
heartbeats. The system is also characterized by employing several digital 
Adaptive Matched Filters (AMFs) on the monitored (heart)Beat Under Test 
(BUT). An AMF is defined by a corresponding (shape) 'template', known as 
'Matching Template', which is mathematically represented by a corresponding 
set of parameters. Different signal patterns/shapes may be identified by AMF 
by correspondingly defining other sets (i.e., templates) of parameters. The 
AMFs are utihzed for allowing the system to detect and classify Viable 
heartbeat' signal, and, after detecting several viable heartbeats (i.e., according 
to a predetermined criteria), to adapt to a person currently being monitored. By 
using the term Viable heartbeat' it is meant to include every possible 
shape/pattern of a signal that resembles common Normal (see Figs, la and lb) 
and common (aberrant) Pathological (see Figs. Id and le) shape/pattern. In 
connection with the EEG signals, a known process and model are utilized. The 
process is the HMM, the parameters of which are continuously updated in order 
to allow the system to converge to the inherent sleep patterns of a patient, the 
sleeping behavior of which is to be analyzed. The model is the R&K model, and 
it provides a chnician a chnical-based set of rules that allow the clinician to 
distinguish between different sleep stages. According to the present invention, 
the clinical-based set of rules have been translated into mathematical process, 
for allowing a fully-automated and accuracy-enhanced analysis of sleep stages 



wo 03/057025 PCT/n.03/00005 
-25- 

and sleeping disorders. In connection with the Respiratory signals, the present 
invention utilizes Fuzzy Logic Decision Process (FLDA), for obtaining more 
reahstic decisions regarding the feasibiUty of Respiratory disorders. 

The present invention is further characterized by utihzing only one (novel) EEG 
electrode, instead of four EEG electrodes, as in conventional systems, said novel 
EEG electrode being attached to a patient in a location that allows, in 
conjunction with the novel processes, to obtain the correct sleep patterns. 

The system is designed to measure biological signals of a patient that could be 
located essentially anywhere, including at his home. Accordingly, data 
representative of the physiological signals of the patient may be sent in real- 
time through a data network, such as the Internet, Local Area Network (LAN) 
or by utilizing a telephone, to a laboratory server, in which the data is processed 
and analyzed, and the resisting data could be transmitted to a Web server, for 
introducing the complete, or selections of, said data, to a clinician according to 
his preferration. 

The disclosed system includes a hardware part, for acquiring the required 
physiological signals, and a software part, for processing and analyzing said 
signals. 

The hardware comprises at least: 

(1) EGG electrode strip, for acquiring data related to the condition of the 
examined heart; 

(2) Forehead EEG electrode stripe, for acquiring data related- to the sleep stages 
of the examined patient; 

(3) Chest and abdomen respiratory belts, for acquiring data related to the 
respiratory status of the examined patient; 

(4) Oximeter, for measuring the Oxygen saturation level in the patient's blood; 
and 

(5) Movement detector being placed on limb of a person, or elsewhere. 
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The software handles tlie following aspects: 
1^ Electrocardiogram (EGG) 

According to the present invention, whenever a person is coimected to the 
system, the system adapts itself to the connected person by performing the 
following steps: 

a) Detecting essentially every signal, that has at least a minimum 
resemblance to a heartbeat-like shape (i.e., a viable heartbeat), that is 
contained within the first monitored ECG segment (i.e., portion), the 
duration of which is predetermined by the clinician, and may be varied by 
him according to the circumstances (e.g., whenever the shape of the 
heartbeats of the monitored person are problematic), by comparing the 
signal to a reference 'Template No. 1', which is initialized with a 
universal reference shape (hereinafter referred to as the 'Universal 
HeartBeat' - UHB, see Fig. la), which was obtained by averaging several 
normal-shaped heartbeats of several persons; 

b) Employing a predetermined criteria for averaging the viable heartbeats 
detected by use of the UHB, for obtaining an average heartbeat shape, or 
appearance, that is unique to the monitored person (hereinafter referred 
to as the 'Patient-Specific Universal HeartBeat' - PSUHB, see Fig. lb). 
Template No. 1 is assigned the characteristics of PSUHB in order to 
utiHze it as the unique reference of the person being monitored; 

c) Comparing additional heartbeat-like signals, in corresponding additional 
portions of the ECG signal, to the PSUHB reference; 

d) Whenever a new viable heartbeat is detected by using the updated 
Template No. 1 (i.e., PSUHB), storing it and analyzing its medical 
implications; 

e) updating Template No. 1 (i.e., PSUHB), according to a predetermined 
criteria, by utilizing additional new viable heartbeats for re-averaging the 
monitored person's unique viable heartbeats, according to predetermined 
criteria, for allowing the system to enhance its adaptation to the 
monitored person; and 
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Viable Heartbeats 

Prior to employing Template No. 1 on a Beat Under Test (BUT), tbe system first 
identifies the relative time-wise location of the BUT within the continuous EGG 

signal. 

The BUTs are identified by manipulating several signals associated with the 
corresponding output channels of the WTA. The number of the channels 
depends on the fi-equency at which the EGG signal was originally sampled. In a 
first WTA output signal (commonly referred to as the 'high scale'), the high 
frequency content of the EGG signal is emphasized (see FIG. 8a). The lower the 
'scale', the lower the dominant frequency of the EGG signal that is emphasized 
(see Fig. 8b and 8 c). 

According to the present invention, the absolute values of the first three 
'scales' (Fig. 8a, 8b and 8c) are 'summed-up', resulting in a signal (Fig. 8d) 
having a high peak, which allows an easy detection of 'R-peak' of every 
heartbeat. Due to the latter maniptdation, the 'R-peak' of a heartbeat wiU 
always be the most perceptible point in the heartbeat, and, therefore, wiU be 
easily detected and will essentially eUminate the risk of misdetection of R- 
peaks (i.e., due to noises, movements, T waves and other types of interference). 

As shown in Fig. 8d, there are two points that have a larger energy, i.e., the 
heartbeats marked as 'AR'. Nevertheless, the lower-energy R- heartbeats are 
viable heartbeats and will be ansJyzed accordingly. From the signal depicted in 
Fig. 8d, the R-peaks are "placed" in the filtered EGG signal, as shown in Fig. 8e. 

The identified R-peaks, as well as the AR-peaks, are utiUzed as timing 
indications, around which the samples (of the EGG signal) are evaluated for 
determining the duration and relative location of each (required) segment (i.e., 
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PQ, QR, RS, QRS and ST). Accordiagly, essentiaUy each BUT is identified by 
the system, after which it is characterized (i.6., its type being determined). 

Each AMF outputs a probability value (between 0.0 to 1.0), indicating the 
resemblance between the BUT to the corresponding template. The larger this 
value, the more the BUT resembles to the corresponding template. The 
corresponding AMF outputs a Normalized Correlation (NC) value, which is 
compared to a first threshold value iVC^j . If iVC > NCjjji , it indicates that a 
viable heartbeat has been detected, in which case it is further processed. 

The PSUHB may be updated, as ECG signal continues to 'feed' the system, 
according to one of the following principles: 

(1) whenever several viable heartbeats are detected within a predetermined 
time interval, the PSUHB is calculated by averaging only the viable heartbeats 
detected in the last ECG segment; or 

(2) as in (1), only that the PSUHB is calculated by averaging the viable 
heartbeats detected in every ECG segment; or 

(3) whenever predetermined number 'n' of viable heartbeats are detected, the 
PSUHB is calculated by averaging only the last 'n' viable heartbeats that were 
detected; or 

(4) as in (3), only that the PSUHB is calculated by averaging gdl of the detected 
viable heartbeats; or 

(5) whenever a new viable heartbeats is detected, the PSUHB is calculated by 
averaging all of the viable heartbeats; or 

(6) whenever a predetermined time-interval elapses. 

Whenever a new PSUHB is calculated, it is used as a new reference template 
for allowing the system to identify new viable heartbeats in future ECG 
segments of same monitored person. As ECG signal continues to be input into 
the system, the system converges from the generalized universal (mean) 
heartbeat patterns into the specific patient's heartbeat patterns (i.e. the system 
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'familiarizes' itself with, or 'learns', the inherent cardiac condition of the 
currently monitored patient). The latter process is employed on Normal 
heartbeat, as well as on Pathological heartbeat, and it allows the TMP to 
determine accurately the cardiac status of the examined person. If the 
examined (current) heartbeat is identified as pathological, a second 
classification process is employed, for determining the exact type of the 
pathology, after which the current identified pathological type is utihzed for 
calculating a new mean pattern that is essentially more accurate than the 
previous calciilated mean pattern. 

At the time a (new) person is connected to the system, the system does not have 
any patient-related data, and therefore it must utilize relatively 'weak' criteria 
to start identifying the first person's BUTs. Therefore, the system utHizes the 
UHB, which is associated with a relatively low threshold value NCtjji. 
However, after the first PSUHB is calculated, i.e., after the first EGG segment 
provides the system with some unique patient-related data (i.e., data that 
reflects the inherent cardiac condition of the person currently being monitored), 
a 'harsher' criteria is employed by the AMF, by utilizing a second threshold 
value NCTH2^i.e., NCiji2> ^Cjjji). 

Whenever a viable heartbeat is detected, which is not pathological heartbeat, 
the system considers it as Normal heartbeat. Each Normal heartbeat is 
segmented into four clinically important segments: QRS, ST, PR and QT. The 
segmentation process is based on the prior knowledge of the relative location of 
the R-wave, T-wave and P-wave of each normal heartbeat, which were found by 
employing the WTA at an earher stage, as described before. 

Accordingly, the system finds the corresponding segments boundaries by: 

(1) Finding the peak of the Q-wave; 

(2) Finding the onset of the Q-wave; 

(3) Finding the peak of the S-wave; 

(4) Finding the end of the S-wave; 
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(5) Finding the peak of the T-wave; 

(6) Finding the onset of the T-wave; 

(7) Finding the end of the T-wave; 

(8) Finding the onset of the P-wave; 

The system automatically searches for the above-specified boundaries wkQe 
synchronized with the R-R interval, which was accurately found by the WTA in 
an earlier stage of the analysis. Accordingly, the system provides a person- 
adaptive and heartbeat rate adaptive segmentation process. 

T-wave detection 

The major problem in the segmentation process is to precisely characterize the 
T-wave, i.e., to precisely detect its shape and phase, which are critical for a 
precise determination of the ST segment and the QT interval. According to the 
invention, the T-wave of a heartbeat is characterized, i.e., a decision is 
automatically made, regarding whether the T-wave is 'normal, 'upside-down' or 
'bipolar'. The T-wave characteristics are found by identifying its onset and end- 
point. 

ST analvsis 

One of the clinically most important parameter associated with cardiac 
condition, is the ST depression/elevation, as this parameter may indicate a 
symptom for heart attack. In order for the system to evaluate the ST segment, 
it selects a portion of the ECG, which is contained within the isoelectric interval 
(the isoelectric inte rval is defined as the time between the end of the former T 
wave and the onset of the current P waveV which is located between two 
consecutive heartbeats. Then, the system 'extracts' 0.1 Second of this interval, 
which is utilized as a reference segment, to which the ST segment is compared. 
An absolute margin is calculated between the mean magnitude of the ST 
segment and the mean magnitude of the reference segment, ST 
depression/elevation is identified if the absolute margin is larger than a 
threshold value. 
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QT analysis 

The duration of the QT interval is calculated, as well as the variability of said 
duration and the corresponding spectrum, by utilizing models known as the 
'AR-models'. 

Sometimes, the 'normal' heartbeat of a person does not include a perceptible R- 
wave (i.e., which the system can not identify or detect). Therefore, the system is 
configured to search for P-waves in these heartbeats, measiire the QRS 
duration and check if the heartbeat is not premature. If these conditions meet 
certain criteria, the system considers the corresponding heartbeats as viable 
heartbeats, and the AMF's template is updated accordingly and automatically 
by the system. 

The system also calculates Heart Rate (HR) and Heart Rate Variability (HRV) 
parameters, which are easily derived from the . ECG signal by employing 
techniques that are well known to skilled in the art. Calculation of HR and 
HRV is implemented by utilizing only Normal heartbeats. 

The system is capable of detecting HR abnormalities from different 
physiological signals such as those derived from the signals associated with a 
heart condition, the patient movement activity, sleep stage, etc., by utilizing an 
automatic physiological events detection and a set of fuzzy logic rules. 

Patholofifical Heartbeats 

Whenever a viable heartbeat is detected (i.e. by use of Template No. 1), the 
system compares it to two commonly (universally) known aberrant heartbeats 
(i.e., PVCl and PVC2), by which Templates No. 2 and 3 are initially assigned: 

(1) Template No. 2 - Universal Pathological HeartBeat No. 1 (UPHBl) 
template. This template is initialized by the shape of the PVCl, which 
represents essentially any signal pattern that resembles a first commonly 
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known aberrant heartbeat (i.e., PVCl, see Fig. Id), and was obtained by 
averaging several aberrant heartbeats from several persons (having the same 
pattern as the first known aberrant pattern); and 

(2) Template No. 3 - Universal Pathological HeartBeat No. 2 (UPHB2) 
template. This template is initialized by the shape of the PVC2, which 
represents essentially any signal pattern that resembles a second commonly 
known aberrant heartbeat (i.e., PVC2, see Fig. le), and was obtained by 
averaging several aberrant heartbeats from several persons (having essentially 
the same pattern as the second known aberrant pattern). 

Employing these two templates results in two Normahzed Correlation values 
(i.e., NCpi and NCp2 , respectively). Normally, NCpi ^ NCp2 . Therefore, the 
system selects the UPHBi (i= 1 or 2) that yields the maximal NC (i.e., between 
NCpi and NCpi), and compares its NC to a third threshold value (i.e., NCjjj^). If 
the condition AC > NCj'jj^ is met, the system decides that the viable heartbeat 
has the corresponding aberrant pattern type (i.e., UPHBI or UPHB2). 
Otherwise, a further analysis is carried out, for identifying the presence of noise 
and/or movement artifacts. The noise and movement artifact detection is based 
on a (specific) patient movement learning process and environmental noise 
elimination process. If no noise or artifacts are detected in the heartbeat, it is 
compared to a forth threshold value NCjjj^, which is selected so that 
NCfu^ < NCffj^ , for essentially eliminating the risk of misdetection of aberrant 
heartbeats. If NC <NCxha, it indicates that the heartbeat is so deformed or 
corrupted that its shape is unknown. 

After detecting several aberrant heartbeats, the system may learn the patient- 
specific aberrant heartbeat patterns and introduce to a cHnician a 

corresponding sequence of pathological heartbeats. 

If the system decides that the viable heartbeat essentially matches the UPHBl 
template, or the UPHB2 template, it stores it for calculating new corresponding 
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templates (i.e., PSUPHBl and PSUPHB2, respectively), which reflect the 
corresponding patient-specific aberrant heartbeat. Calculating PSUPHBl and 
PSUPHB2 is implemented according to the same principles described in 
connection with the PSUHB. 

Pathological Sequences 

After a viable heartbeat is identified as aberrant heartbeat, it is stored, and a 
corresponding sequence of pathological heartbeats is established therefrom. 
There are several commonly known pathological sequences ( 'A' and 'N" indicate 
'Aberrant' and 'Normal' heartbeats, respectively): 

(1) VT - A, A, A, A, A, (generally seen in a tachycardia-type disorder 
associated with the ventricles). VT is detected by finding a sequence of 
pathological heartbeats and instantaneous Heartbeat Rate (HR) that exceeds a 
clinically predetermined threshold value; 

(2) SVT - PAC, PAC, PAC, PAC (generaUy seen in another type of 
tachycardia disorder); (PAC - Premature Atrial Contraction, as seen in the ECG 
as a normal heartbeat morphology only if it is premature); 

(3) Bigeniny - A, N, A, N, A, N, ..,; 

(4) Trigeminy 1 - N, N, A, N, N, A, N, N, A, 

(5) Trigeminy 2 - A, A, N, A, A, N, A, A, N, 

(6) Tachycardia - accelerated heartbeat rate. Tachycardia is detected by finding 
several consecutive heartbeats having a rate that is above a certain clinically 
predetermined threshold value; and 

(7) Bradycardia - too low heartbeat rate. Bradycardia is detected by finding 
several consecutive heartbeats having a rate that is below a certain clinically 
predetermined threshold value. 

2) Automatic Sleep Staging (ASS) system 

The ASS system is configured to identify the following sleep stages: 
wakefulness, sleep Stage 1, 2, 3 and 4, and Rapid Eye Movement (REM), which 
is associated with Dreaming stage. The ASS system is also capable of 
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identifying other special events, such as K-complexes, spindles, alpha 
intrusions, body and electrode movements, eye blinks, awakenings, arousals 
and micro-arousals by utilizing only one electrode. 

The ASS system performs, in general, two analyzing processes. The first process 
is associated with 'on-line' (i.e., 'real-time') process. Execution of this process 
involves real-time extraction of several features and evaluation of the features 
every predetermined time interval (e.g., every 30 seconds), in order to reach a 
decision regarding the current sleep stage. The second analyzing process, being 
an 'off-Une' process, is employed (normally) after several hours of sleeping, for 
allowing the ASS system to converge to (i.e., 'learn') the inherent sleep behavior 
of the monitored person. 

The first (i.e., 'on-line') analyzing processes involves the following steps: 

(1) Filtering the EEG signal, EOG and the EMG signals; 

(2) Adaptively segmenting the resulting filtered EEG and EOG signals into 
QSSs; 

(3) Extracting the required features from the resultiug quasi-stationary EEG, 

EOG and EMG segments; 

(4) Preliminary classification of each EEG QSS. 

(5) Whenever relevant, detecting special temporal events contained within the 
EEG QSS; 

(6) Special events classification, by using environment characteristics; 

(7) If an EOG signal exists, detection of EOG REM, by employing related 
process; 

(8) Real-Time determination of sleep stages; 

(9) Real-Time determination of awakening, arousals and micro-arousals; and 

(10) Repeating steps (1) to (9). 

The second analyzing process involves the following steps: 



(11) Initial division of the firequency groups into sub-groups; 



wo 03/057025 PCT/n.03/00005 
-35- 

(12) Employment of Baum- Welch on the HMM model; 

(13) Refining the special event classification of step (6) in the first analyzing 
processes. 

(14) EMG energy 'lows' (i.e., minimal points) detection process using 
hierarchical clustering process that is based on fuzzy logic; 

(15) Refining sleep staging process; and 

(16) Refijiing awakening, arousals and micro arousals detection. 

(1) Filtering 

The raw EEG signal is filtered by two types of (digital) Band-Pass Filter (BPF), 
which are implemented by DSP technique commonly known as 'Finite Impulse 
Response' (FIR). The first FIR-based BPF (l-35Hz) cancels drifts in the basehne 
of the raw EEG signal and acquires the desired EEG bandwidth. The resulting 
signal is hereinafter referred to as the 'filtered EEG'). The second FIR-based 
BPF (20-35Hz) acquires the EMG signal (the tonus muscle is very close to the 
Fp electrode). The filtered EEG is decomposed into two signals by utilizing a 
third and a fourth FIR-based EPFs. The third and fourth BPFs are referred to 
as the 'Slow Frequency EEG' (SFEEG) and the 'Fast Frequency EEG' 
(FFEEG), respectively. If the EGG and EMG signals also exist, the 
corresponding drift in the EEG signal is canceled by employing corresponding 
BPF(s). 

(2) Adaptive segmentation 

EEG signal is a non-stationary process, namely, its statistical characteristics, 
such as frequency and moments, are time-dependent. In order to obtain a 
reliable identification of sleep stages of a monitored person, a first step is 
implemented, in which the EEG signal is partitioned into QSSs using an 
adaptive segmentation process, hereinafter referred to as the 'Adaptive 
Segmentation Process for EEG signals' (ASAEEG). 

The ASAEEG process utilizes the known (i.e., to those skilled in the art) 
Generalized Likelihood Ratio (GLR) measure, for distinguishing between 
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statistically QSSs. For example, given two time series, the GLR meastire 
compares the statistical properties of the two series and outputs a value, which 
represents the statistical difference (or 'distance') between the two series. For 
example, a zero value means that the two series have exactly the same 
statistical properties, a fact which leads to the conclusion that the two series 
belong to the same statistical process. 

The ASAEEG process is characterized by its capability to dynamically change a 
'time-window' (hereinafter referred to as the 'Growing Reference Window' - 
GRW), within which a corresponding portion of the EEG signal is analyzed. A 
signal, which is contained within said portion, usually comprises a dominant 
wave (i.e., which characterizes said portion) and several frequency components 
that are associated with other waves and superimposed on said dominant wave. 

The GRW changes automatically in duration, as well as in its relative location, 
untU an optimized decision with respect to specific examined EEG portion, is 
obtained by utilizing the GLR measure. 

The probabiUty values of the five waves form a vector, which characterizes the 
ciirrent (stationary) segment being analyzed. As additional consecutive 
segments undergo the same process, additional corresponding vectors are 
formed. Every vector is included in a matrix, which characterizes the sleep 
status of the monitored person. The system may automatically update the 
matrix (i.e., by adding additional vectors) according to some criteria. For 
example, the matrix may be updated per time-interval. The larger the matrix, 
the more the system is adapted to the inherent patient's specific sleep behavior, 
resulting in more enhanced accuracy of the corresponding decisions made by the 
system 
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Identifidng features in a quasi-stationary seement 

An EEG signal normally comprises five major types of signals (generally 
referred to as waves): (1) ^-wave (1.0-3.5 Hz), (2) ,9 -wave (4.0-7.0 Hz), (3) a- 
wave (7.5-12 Hz), (4) c-wave (12-15 Hz); and (5) j5-wave (15-35 Hz). 



For each QSS that is identified, several, or all, of the following features 
associated with are obtained: 

1) Relative energy in the S -wave Band; 

2) Relative energy in the & -wave Band; 

3) Relative energy in the a -wave Band; 

4) Relative energy in the a -wave Band; and 

5) Relative energy in the -wave Band; 

6) Statistical data, such as 'Mean value', 'Variance' and 'Skewness'; 

7) Duration of the QSS; 

8) Fimdamental frequency of the QSS. 

Data related to features Nos. 1) to 8) are obtained in connection with the 
filtered EEG, as well as the EMG signals, while features Nos. 1) and 2) are 
obtained in connection with the EOG signal. 

(4) Preliminary classification of each EEG quasi-stationary segment . 
After a QSS is characterized (i.e., according to step 3), it is associated with one 
of the three fi-equency groups (i.e., HFG, MFG or LFG), thereby creating three 
corresponding clusters, by employing fuzzy kgic rules on the extracted features. 

A real-time decision, regarding a current sleep stage, is based on evaluating 
several consecutive QSSs. Therefore, the system automatically makes sleep 
stage evaluations after each predetermined time-interval, in which several 
QSSs are contained. For example, a time interval of 30 seconds has been found 
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to be wide enough for containing several QSSs, wMle still sustaining the 'real- 
time' attribute of the system. 

(5) Detecting special temporal events contained within the EEG quasi- 
stationary segment. 

The special events 'K'-complexes, body and electrode movements, rapid eye 
movements and eye blinks are temporal events, and are detected by exploiting 
the fact that they differ from the ongoing EEG signal by their energy, skewness, 
duration and fundamental frequency. 

(6) Special events classification. 

Because there is a strong correlation between special events and their adjacent 
signals, an 'environment enrichment' principle is employed, namely, each 
special event is classified after characterizing the frequency content of the 
corresponding adjacent QSSs, and based on its temporal morphology. For 
example, if a special event is detected, which has a distinguishable temporal 
morphology, and its adjacent QSSs belong to the MFG, it is most probable that 
this event is a 'K-complex'. 

After the time-wise (temporal) events are identified, they are essentially 
removed from the ongoing EEG signal for allowing the system to analyze its 
corresponding spectral features (i.e., the relative energies of the S , &, a and 
the J3 waves) by utilizing the HMM models. 

n") Detection of EOG REM 

REM signals are associated with dreaming stages, and are detected by 
exploiting the fact that REM are paradoxical and have higher energy regarding 
their environment. 

(8) Real-time sleep staging 



Detecting changes from one frequency group to another 
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The ongoing EEG signal is substantially a Markovian process, namely; the 
current status of a sleeping person largely depends on its past behavior, a fact 
which allows employing the Hidden Model (HMM) statistics model on the 
ongoing EEG signal. 

In order to simplify the problem of identifying changes between the three 
frequency groups (i.e., HFG, MFG and LFG) a 'Three-States' model (TSM) is 
utOized, which is based on the assumption that there are three major frequency 
activities (or 'states') involved in a sleeping person. There are various 
probability values associated with passing from one frequency group (e.g., HFG) 
to another frequency group (e.g., LFA). The task of the system is to reach a 
decision regarding the current stage, by evaluating the (known) predecessor 
stages and by utilizing statistics-based estimation process. 

Accordingly, the above-mentioned five sleep waves have been grouped into three 
frequency groups (stages). Each frequency group corresponds to a different 
sleep stage (i.e. hght sleep, dreaming and awakening). The decision-making 
process, regarding the sleep stages, is based on the HMM model. 

Each stage is characterized by a 4-dimensional (probabiHtjO gaussian. 
Accordingly, a two-dimensional probability graph is formed, in which the 'X' 
axis is the probability that the current segment is the a -wave, and the 'Y' axis 
is the probability that current segment is the S -wave. The system initializes 
the probability graph with three 'global' mean 'staxting-points', each of which 
represents a different frequency group. The global starting-points are obtained 
by averaging EEG starting points of several persons. Alternatively, individual 
global starting-points are calculated, by employing a clustering process, which 
is based on fuzzy logic, on the pre-determined frequency groups. 

Sleep stasdng 

The final decision, regarding the relative time-wise location and duration of the 
sleep stages, is obtained by implementing the following steps: 
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1) Averaging each one of the three clusters, i.e., the QSSs that were found in 
the last 30-seconds time-interval, that are associated with the same frequency- 
group; and 

2) Employing a set of R&K rules on the averaged clusters and on the temporal 
events. 

(9) Awakenings, arousals and micro-arousals detection: 

A fuzzy logic process, based on the rules for arousal detection, is used, in order 
to automatically detect awakenings, arousals and micro-arousals. 

(10) Repeating steps (1) to (9) for each new 30-seconds time-interval . 

(11) Initial division of the frequency groups into subgroups. 

Since each person has essentially different sleep patterns, which are uniquely 
represented by corresponding signals, it is allowed to characterize the person's 
sleep behavior more accurately, by dividing the frequency groups into 
subgroups according to unique physiological properties/criteria that 
characterize the corresponding frequency group. For example, the HFG is 
divided into two subgroups - a fii-st subgroup is characterized by having 
increased (^-wave activity, and a second subgroup is characterized by having 
increased ^ -wave activity. Likewise, the MFG is divided into three subgroups. 
The LFG remains as one group. 

(12) Employment of Baum -Welch on the HMM model. 

A Gaussian ProbabiHty Function (GPF) is employed in connection with each of 
the (frequency) subgroups. Each 'gaussian' is characterized by two parameters, 
namely, its 'center' (i.e., the 1.0 probabUity) point and variance, and represents 
the probability of a given quasi-stationary segment to belong to the 
corresponding subgroup. 

The Baum- Welch process is implemented in the second analyzing process (i.e., 
in the 'off-line' process), and its task is to update the gaussian functions. The 
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Baum-Welch process is employed every several hovirs (i.e., of sleeping) on the 
last known GPFs and the HMM model, in order to update the GPFs. The more 
times the Baum-Welch process is employed, the more accurately the ASS 
system converges to the inherent (unique) sleep patterns of the monitored 
person. 

(13) Refining special events classification. 

Since the 'real-time' decisions are hased on a global GPF (i.e., without 
considering 'learning' process), these decisions might be inaccurate. However, 
the updated GPFs (i.e., person's specific GPFs) allow the ASS to reprocess (i.e., 
'refine') the data on which the 'real-time' decisions were based, and obtain 
accurate decisions regarding special events. 

In other words, the ASS system adapts itself to the actual three groups of 
frequencies that xmiquely represent the inherent sleeping patterns of the 
monitored person. Each new detected quasi-stationary segment of the EEG 
signal contributes to the decision making process, as the accumulated effect of 
all the quasi-stationary segments allows the system a better 'learning' of the 
patient. 

(14) EMG energy 'lows' (i.e.. minimal energy points) detection process . 

EMG energy 'lows' are detected by using hierarchical fuz2y clustering process 
which recursively partitions the low EMG energy group into a low energy group 
and a high one, until stopping condition is fulfilled. 

(15) Refining sleep stages classification . 

Refining sleep stages classification is carried out essentially the same way as 
the special events classification (see step 13). 

(16) Refining arousal detection . 

Arousal detection is carried out essentially the same way as the special events 
classification (see step 13). 
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3) Respiratory 

A probability-based decision, associating Respiratory disorder with an apnea (or 
hypoapnea), is made by the system after it considers data from five data 
sources/channels: (1) EMG, (2) Limb movement, (3) Chest effort, (4) Abdominal 
effort and (5) Blood saturation. In order to simplify the decision-making process, 
the five data sources/channels have been grouped into three groups: (1) Effort, 
which includes the Chest Effort signal and the Abdominal Effort signal, (2) 
Movement, which includes the EMG signal and the limb Movement signal 
and (3) Oxygen Saturation, which is measure by utilizing an Oximeter. 
Corresponding processes have been developed for extracting relevant features 
from each of the above-described groups. 

A normal breathing signal is characterized by having a low frequency signal, 
the fe-equency and magnitude of which are essentially constant. In this case, the 
envelope of this signal does not change with respect to time. An hypoapnea is 
clinically defined as having the respiratory airflow decreasing by at least 20% 
fi:om its initial value, for duration of at least 10 Seconds (Amy L. Meoli et al. 
"Hypopnea in Sleep-Disordered Breathing in Adults" Sleep, vol 24 No. 4, 
2001). Accordingly, in order to decide whether a Respiratory disorder is an 
hypoapnea (or an apnea, which is clinically more severe, but easier to identify), 
the envelope of the effort signal is first analyzed, by carrying out the following 
steps: 

(1) Detection and Smoothing of the Envelope (DSE) of the effort signal. In order 
to produce the required envelope, the amphtude of each breathing cycle is 
identified, and the shape of the envelope of the breathing signal is 'smoothed' by 
utiUzing a set of fuzzy logic rules, in order to remove irrelevant maximum 
points. 

(2) Peak Detection and Maximum Setting (PDMS) of the envelope signal. After 
establishing the corresponding envelope signal, its maximal points are identified 
and utiHzed as time-wise boundaries for segmenting the effort signal into 
corresponding segments. If there is a group of at least two maximal points (i.e., 
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in the envelope), which are too close (i.e., with respect to time) to one another, 
only the larger maximal point in that group is selected as time houndary. The 
identified time boundaries are also utilized for segmenting the envelope into 
segments that overlap the corresponding effort signal segments. 
(3) Segment Evaluation Process (SEP). Each segment (i.e., in the effort signal) 
is evaluated for determining whether there is more than 20% decrease in the 
magnitude of the effort signal contained within said segment. If there is a 20% 
decrease in said magnitude (i.e., for a duration larger than 10 seconds), a Fuzzy 
Logic Process (FLP) is employed, for evaluating other relevant events, such as 
Oxygen de-saturation, arousal(s), movement(s). 

However, if a decrease less than 20% is not evident in said magnitude, it is 
necessary to determine the 'direction' of the envelope of the evaluated segment, 
i.e., whether it tends Upwards (i.e., its time derivative is positive in each point 
of the corresponding envelope), Downwards (i.e., the envelope has negative time- 
derivative) or it maintains essentially the same magnitude. In order to 
determine the direction of each envelope's segment, a mathematical 
approximation to the envelope's segment is made, by utiHzing a polynomial of 
second order (parabola). Then, mathematical manipulations are carried out for 
identifying the required corresponding time-derivatives of the corresponding 
parabola. Each envelope segment, which does not meet the '20% magnitude 
decrease' criteria, is assigned a corresponding parabola. 

A Segment Combining Process (SCA) is employed on segments that do not 
comply with the 20% of the EPA criteria (see par. 2), in order to adjoin two 
adjacent segments to one (grand) segment. The SCA could be the Viterbi Process 
- which wiU not be discussed, since it is known to those skilled in the art 

(Reference could be made to the web site - http:/www.comp.leeds.ac.uk/scs- 
only/teaching- materials/ http:// www.comp.leeds.ac.uk/scs-only/teaching- 
materials HiddenMarkovModels/ html dev/ viterbi nrocess/sl pg3.html) . 
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(4) Decision Making Process (DMA). According to the present invention, the 
Respiratory system is characterized, by emplojring Fuzzy Logic rules for allowing 
the system to reach more realistic decisions regarding the nature of the 
Respiratory disorders. Accordingly, the DMA is based on a Fuzzy Logic Process 
(FLP), which is utilized for estimating the probability of a Respiratory disorder 
being an apnea, or hypoapnea, while taking into account data from the above- 
described data sources/channels. The FLA gives mathematical interpretation to 
known clinical set of rules (accordiag to which apnea is detected), and, thereby, 
allows the system to mimic the analysis process that is currently made 
manually by a clinician. 

The FLA utilizes several trapezoid-shaped 'ProbabiUty-Wise Windows' (PWWs). 
Each PWW is associated with a corresponding event that is encountered in the 
corresponding data source selected from the group of data sources (1), (3), (4) 
and (5), which are described hereinabove. The term 'event' refers to noticeable 
change that is detected in the corresponding signal, such as sudden movements 
and/or exerted efforts. For example, in connection with the first data source, i.e., 
"(1) EMG", an event means that a rise in the EMG level has probably been 
detected. Each PWW has a duration that is relatively much larger (i.e. wider) 
than the duration of the corresponding event. In addition, each PWW comprises 
a positive slope portion (representing a gradual increase from 0.0 to 1.0 
probabihty value), a 1.0 probabUity value portion and a negative slope portion 
(representing a gradual decrease fi:om 1.0 back to 0.0 probabihty value). 

The time duration of each part in each PWW has been optimized according to 
medical criteria, as well as its relative location (i.e., time-wise synchronization) 
of each PWW. For example, if the monitored person starts suffocating (i.e., the 
magnitude of the effort signal decreases by more than 20%), a corresponding 
delayed decrease in the Oxygen saturation level in the person's blood is 
expected to be observed. A normal delay (i.e., between suffocation and noticeable 
decrease in the Oxygen saturation level) is of a couple of seconds (e.g., 4-6 
seconds). Accordingly, a first step wiU take place, according to which an Oxygen- 
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related PWW will 'open' (i.e., placed) in a way that its (vertical) center is 
aligned with the ending point of the apnea event. Then, a second step will take 
place, according to which the system analyzes the time location of the Oxygen 
event in relation to the Oxygen-related window. At the end of the latter 
analysis, the oxygen is assigned a probability value. For example, in an ideal 
apnea event, the delay between the Oxygen event and the suspected apnea 
event matches exactly the cUnically predetermined delay, in which case the 
instant of the Oxygen event matches exactly the time center of the Oxygen- 
related PWW, in which case the Oxygen event is assigned the probabihty value 
1.0, Accordingly, the probabihty that the suspected apnea is a 'real' apnea is 
1.0. However, the delay may vary from one person to another and from one 
apnea to another, in which case the corresponding Oxygen event will be 
assigned a probability value other than 1.0. In some cases, the delay will be 
such, that the Oxygen event will be timely located outside the Oxygen-related 
PWW, in which case it will be assigned the probability value 0.0. Accordingly, 
the suspected apnea will not be regarded as a 'real' apnea (i.e., the essential 
condition of a following decrease in the Oxygen saturation level - in the 
expected time window and time location - is not met). 

Phase Detection Process (PDA) 

A change in the relative phase between the two effort signals (i.e., the chest 
signal and the abdominal signals) may indicate an obstruction in the upper 
airways of the monitored person. The PDA performs the following steps: 

a) Generating a signal that represents the phase between the two effort signals. 

b) Whenever applicable, finding minima and maxima, the duration of which is 
more than 10 Seconds. 

c) Whenever applicable, finding adjoining events. 
Movement Detecting Process HVIDP) 

A movement is characterized by a noticeable increase in the signal energy. The 
task of the MDP is to allow distinguishing relatively large movements from the 
surrounding environment interferences. The MDP performs the following steps: 
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a) CEilculating the signal energy in overlapping time 'windows' with length L. 

b) In each segment, finding the points with maximum values. 

c) Calculating the normahzed distance between each two of the points found in 

step b). 

d) Finding the Maximal Normalized Distance (MND) in the normalized 
distances calculated in step c). 

e) Determining whether there is a movement, by comparing adjacent MNDs as 
found in step d) 

Oxygen Saturation Event Process (BSEA') 

A blood saturation is defined as an 'event' whenever there is at least 4% 
decrease from the last steady level. The BSEA performs the following steps: 

a) Smoothing the signal with an averaging window of length L. 

b) Finding local minima points in the smoothed signal. 

c) Comparing the level of the minimtun point to the level of the minimum's 
base. 

d) Classifying events into one of two categories, according to a threshold level 
being lower than 4%: (1) events in which a magnitude decrease has been 
detected that is less than the threshold, and (2) events in which a magnitude 
decrease has been detected, which is more than the threshold. Events associated 
with the second category are assigned a higher weight in the final decision 
making process than the first category's events. 

The present invention is directed to a method for measuring a pluraUty of 
physiological signals and output a data, being essentially noiseless and artifact- 
free, the content of which is correlated with corresponding physiological data. 
The present invention is also directed to a system capable of identifying 
'normal' physiological activities of individual patients and adjusting itself 
according to the identified normal physiological activities. The operation 
principles of the system may be better understood with reference to the 
drawings and accompanjdng description. 
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Fig. 2 schematically illustrates the initial identificatioii of signals having 
essentially a heartbeat-like waveform, according to a preferred embodiment of 
the present invention. In order to avoid wasting processing time on signals that 
do not resemble heartbeats, a first classification process is employed (202) on 
the current Beat Under Test (BUT) 201. 

Even though there is a considerable variance involved in sign^s that represent 
heartbeats of different persons, these signals are characterized by having 
typical (noticeable) patterns. An average heartbeat pattern was estabhshed by 
averaging heartbeats of several persons (the 'average heartbeat pattern' is 
hereinafter referred to as Universal HeartBeat - UHB). Accordingly, the first 
monitored BUT is compared to the UHB (203) in order to decide whether it is to 
be considered as a viable heartbeat. The comparison process and resulting 
decision (200a) is implemented by utilizing a digital filter that is commonly 
referred to as a "Matched FQter", which allows implementation of a "Matched 
Template"). The Matched Miter (MF) is configured to identify predefined signal 
pattern (i.e., "template") by employing a corresponding set of parameters. The 
resxdting decision 203 of the MF 200a involves a correlation factor, the value of 
which may be between 0.0 and 1.0. The higher the value of the correlation 
factor, the more the examined pattern resembles the predefined template of the 
MF. Since the magnitude of the BUTs may have a significant variance (i.e., due 
to movements of the sensing electrode, clinical reasons and non-standard 
amplifiers/equipment), a phenomena that may iU affect the resulting decision, 
the BUT is Normahzed (not shown). MF 200a allows identifying even severely 
deformed heartbeats (i.e., pathological heartbeats). However, if the current BUT 
is too corrupted/deformed (i.e., there is very low correlation between the BUT 
and the MFs template), it is not processed for detecting pathologies. 
Nevertheless, the corrupted heartbeat may be stored for (whenever required) 
allowing performing future analysis. 

Assuming the system identifies heartbeats that essentially match the UHB (i.e., 
'proper heartbeat'), it stores these heartbeats, and, according to a first 
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approach, whenever a predetermined time interval elapses (204), calculates an 
average pattern (205). The averaged pattern is the 'Patient-Specific UHB' 
(PSUHB) 206, which is stored in memory 208 (or in memory 207). If the BUT 
does not match the UHB (203), the system checks the next BUT (201). 
According to a second approach, the PSUHB is recalculated each time the 
system identifies new 'n' proper heartbeat (204). Every time a new PSUHB is 
calculated, it is utiHzed as a new universal heartbeat reference (203), to which 
the new BUTs are compared. The PSUHB is constantly updated, providing that 
new BUTs are sampled, thereby allowing the system to adapt itself to (i.e., 
learns) the patient's inherent cardiac condition. The learning process allows the 
system to reach accurate decisions regarding the patient's cardiac condition. 

Each time a proper heartbeat is identified, it is stored (207) and analyzed for 
further classification, e.g., for identifying pathological type (209), if such exists. 

Eig. 3 schematically illustrates the initial identification of pathological-like 
waveforms, according to a preferred embodiment of the present invention. After 
a BUT is identified as UHB (203), it is compared to Universal Pathological 
HeartBeat (UPHB) 302. The comparison process in this case is implemented 
also by utilizing a MF, which has been configxired to identify particidar 
pathological patterns. 

If a current BUT meets the criteria of UPHB (302), it is further processed in 
order to determine its exact pathology type (305). In addition, a Patient-Specific 
Universal Pathological HeartBeat (PSUPHB) (307) may be determined in a way 
similar to determining PSUHB (206), i.e., steps 303 and 306 (Fig. 3) are 
essentially identical to steps 204 and 205 (Fig. 2), respectively. Each time a 
PSUPHB (307) is calculated, it is utilized as a new reference pattern (i.e., 
template) for Matched Filter 302, to which new UHB (203) are compared, until 
a new PSUPHB is calculated again (307). This way, the system adapt to (i.e., 
'learns') the patent's inherent pathological cardiac condition. 
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Current UPHB (302) and PSUPHB (307) may be stored in memory 304 and/or 
in memory 308. 

Fig. 4 schematically illustrates the adaptation principle of the system to the 
Patient-Specific Universal HeartBeat (PSUHB), according to a preferred 
embodiment of the present invention. The adaptation process continues as long 
as new viable heartbeats are identified as having high normalized correlation 
(NC) value, as described hereinafter. 

After detecting the j&rst R-R interval (404), the first BUT is isolated (401). In 
the first iteration the system correlates the first BUT with the Universal 
HeartBeat (UHB) stored in the system and a 'Normalized Correlation' (NC) 
value is calculated (402), i.e., by utiUzing a Matched Filter (not shown) in the 
way described hereinabove. Viable heartbeats having high correlation values 
are averaged, whether on a time interval basis or viable heartbeats count, and 
the new PSUHB (405) is utilized as a reference (i.e., MF template) for the next 
heartbeat(s). Namely, if the condition NC > NCj-ffi (403) is met, wherein NCjjji 
is a threshold NC (having a typical value of 0.7), a new PSUHB is calculated 
(405), and the new PSUHB is stored in the storage array (406a) as the new 
reference pattern/template. Along with storing the (new) PSUHB, the next BUT 
is isolated (405a) and the process repeats itself. The more viable heartbeats are 
identified (i.e,, having high correlation values), the more accurately the system 
converges (i.e., adapts) to the specific patient being examined. However, if a 
second situation occurres, i.e., a patient normally has heartbeats that do not 
conform to the reference pattern/template stored in the storage array (i.e. the 
condition indicated 403 is not met for relatively large duration), the last-known 
viable heartbeat, having the maximum correlation value (LB{NC) \^^ ), wUl be 
regarded as the current PSUHB (404). This way, the system is capable of 
adapting itself to a patient even in cases in which the patient has abnormal (but 
non-pathological) viable heartbeat patterns, even though such cases are rare. 
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If a BUT does not have a P-wave, it is likely to be regarded as aberrant 
heartbeat. Otherwise, the BUT might be a special 'normal' heartbeat, which 
characterizes the heart activity of specific person. In the latter case, it is risk 
that the BUT wiU not be identified as a viable heartbeat. In order not to allow 
midsection of such BUTs, the P-wave of the BUTs is searched for (407) and if 
found, the patient's problematic BUT would be considered as his unique 'viable' 
heartbeats. If a P-wave is found in the current BUT, the current BUT is utilized 
for calculating the current PSUHB (405), otherwise, the BUT is regarded as 
severe pathological heartbeat, or it assumed that the patient is having a critical 
cardiac condition (408). 

In a forth situation (not shown), the heartbeats are very corrupted for a 
relatively long period. In this case, the system will halt, asking the cUnician 
whether to continue the analysis after selecting heartbeat patterns of the 
patient as viable heartbeats, or to abort the analysis. 

Fig. 5 schematically illustrates the heartbeat classification process, according to 
the preferred embodiment of the invention. The function of blocks indicated 500 
and 501 is similar to these of blocks indicated 400 and 401 described in 
connection with Fig. 4. The PSUHB, which is calculated per patient according to 
the process described in connection with Fig. 4, is utihzed as a reference 
template. The current BUT is correlated with the PSUHB (502) that is stored in 
a storage array (not shown), and a NC is calculated (502). 1£ the condition 
NC>NCxH2 is (502), wherein NCjjj2 is a second threshold of NC, the 
current BUT is considered to be a viable heartbeat. Since at this stage the 
system is already adapted, at least to some extent, to the patient being 
monitored (i.e., it utilizes the PSUHB rather than the UHB), the criteria, which 
is utilized by the system for identifying viable heartbeats may be harsher, 
namely, NC^hi > NCjjji . The latter threshold reflects the large variance, and 
therefore low correlation value (e.g. 0.7), of heartbeats of different people, and it 
is required only at the first stage of the analysis, in which the system adapts 
(i.e. converges) itself to specific patient. After completing the adaptation phase, 
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and due to the relatively small variance (and, therefore, high expected 
correlation value) of the heartbeats of the same specific patient, the value of the 
NC could be increased (e.g. to 0.85), for ensuring high precision based decisions 
regarding the classification of BUTs. 

If the condition 503 is not met, i.e., the BUT has small correlation value, the 
BUT is likely to contain artifacts interference or noise. Therefore, a second 
decision is made in accordance with movement data (504). If no movement has 
been detected while sampling the cvixrent BUT, this BUT is suspected as 
pathological heartbeat (506), and the suspected heartbeat undergoes another 
process (not shown), the function of which is to classify the exact pathology type 
associated with this BUT. However, if an artifact has been identified, a decision 
is made (504) that the BUT was sampled vmder very noisy environment 
conditions, and it is to be discarded (507). 

Fig. 6 shows testing restdts of the first heartbeat classification process 
according to which viable heartbeats are classified, according to a preferred 
embodiment of the present invention. The system is capable of distinguishing 
noisy viable heartbeat firom noiseless viable heartbeats. Channel (a) shows the 
ECG signal, in which the viable heartbeats have been identified and marked by 
the system (i.e., by corresponding vertical arrows). These viable heartbeats were 
originally marked by yellow vertical lines (on the computer's screen) and have 
been replaced by said vertical arrows, for demonstration purpose. In addition, 
the problematic ECG signal in sections Tl and T2 was originally marked by a 
conspicuous color in order to attract the clinician's attention. Channel (b) shows 
the ECG signal, in which the viable noisy heartbeats have been identified and 
marked by the system (i.e., by corresponding vertical arrows), and channel (c) 
shows the movement that was detected in time intervals T3 and T4, which 
essentially overlap time intervals Tl and T2, respectively. As can be seen in 
Fig. 6, the system is capable of identifying viable heartbeats even in cases in 
which there are relatively large movements of the patient and/or the sensing 
electrode. 
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Fig. 7 schematically illustrates the pathological heartbeat first classification 
process, according to the preferred embodiment of the invention. Generally, 
there are two typical types of aberrant heartbeat patterns associated with 
common heart pathologies. These types of aberrant heartbeat patterns are 
represented by the Universal Pathological heartBeat No. 1 (UPBl) and the 
Universal Pathological heartBeat No. 2 (UPB2), and utilized by the system as 
reference templates. After deciding that a BUT is suspected as pathological 
(701) (see reference numeral 305 in Fig. 3), the system further examines it by 
comparing it to the UPHBl (702) and to UPHB2 (703), and a first and a second 
Normahzed Correlation (NC) values, i.e., NCpy and NCp2, respectively, are 
calculated (702 and 703, respectively). The maximal NC value is chosen (704), 
which associates the BUT with the most probable corresponding pathological 
type of the BUT. If the condition NC>NCth^ is met (705), it indicates that the 
ctmrent BUT is indeed pathological, the type of which is determined according 
to step 704, and the system utilizes this BUT for learning the patient's 
individual heartbeats, and for analyzing possible pathological sequences of 
heartbeats (707). However, if the latter condition is not met, the noise content of 
the current BUT is analyzed, and a decision is made by the system, whether the 
current BUT is noise (709) or another (weaker) criteria should be employed (i.e., 
in order to avoid mis-detection of severely corrupted pathological heartbeats). 
Accordingly, if the condition NC>NCj-ff4 is met (710, NCjjj/^f^ 0.25), the current 
BUT is considered as pathological (706). Otherwise, the system decides that the 
pattern of the BUT is tmknown (711). 

After a BUT is pathologically classified (706), it is utiKzed for updating the 
corresponding Universal Pathological heartBeat (UPB) reference template, to 
which the new BUTs are compared, by averaging the 'n' last known pathological 
heartbeats (of same type). The updated UPB reference template becomes the 
Patient-Specific Pathological Heartbeat (PSPHB) reference, since the system 
converges into (i.e. learns) the pathological patterns of the patient under 
examination. 
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Fig. 8 illustrates R-peak detection by employing the Wavelet process on EGG 
signal, according to a preferred embodiment of the present invention. 
Fig. 9 schematicaUy illustrates the segmentation process of Normal heartbeats, 
according to a preferred embodiment of the present invention. The Wavelet 
Transfer Process (WTA) is utilized for precise identification of the R-peak of 
each BUT. The WTA is also utilized for precise identification of the P-peak, T- 
peak, T(onset) and T(off) of each BUT. The system utilizes the R, P and T peaks 
for identifying also the P(onset), Q-peak, Q(onset), S-peak and S(off) of the BUT. 
In steps 1, 2 and 3, the system analyzes the section on the left-hand side of the 
(known) R-peak, for identifying the Q-peak, Q(onset) and P(onset). In steps 4, 5, 
6 and 7, the system analyzes the sections on the right-hand side of the R-peak, 
for identifying the S-peak, S(off), T(onset) and T(end). Accordingly, the system 
identifies the QRS complex, ST-segment, QT-segment and the PR-segment of 
each BUT, provided that the BUT is identified as a viable heartbeat (as opposed 
to a BUT that is too deformed/distorted or noisy), which allow identification of 
heartbeat abnormalities other than the two commonly known Pathological 
heartbeat shapes. 

Since the system's search for the above-described P, Q, R, S and T points is 
synchronized to the R-R interval (i.e., the duration between two consecutive R- 
peaks), it is clearly obvious that the segmentation process adapts itseK to the 
heart-rate of the monitored patient. 

Figs. 10a to lOe show testing results that emphasize the various segments of 
Normal heartbeats, according to the preferred embodiment of the present 
invention. In Pig. 10a, the 'ST segments are distinguished fi:om the other 
segments by corresponding horizontal lines. Likewise, in Fig. 10b, the 'QRS' 
segments are distinguished from the other segments by corresponding 
horizontal fines, and in Fig. 10c, the 'QT' segments are distinguished fi'om the 
other segments by corresponding horizontal fines. The latter segments were 
originally marked by red color, which is could not be visible in the figures. 
Therefore, the red color has been replaced by horizontal lines to indicate the 
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corresponding segments. Fig. lOd shows the 'QTV segments, and Fig. lOe 
shows the spectrum of the QTV segment. Pathological heartbeats are not 
segmented, as is shown in Pigs. 10a, 10b and 10c. As can be seen in Figs. 10a, 
10b and 10c, the pathological heartbeats is not segmented. 

Figs. Ha to llh show testing resxilts that were obtained by monitoring four 
persons having different types of pathological heartbeats, according to the 
preferred embodiment of the present invention. Fig. 11a shows 'VT type 
pathological heartbeats that were taken from a first patient. Fig. lib shows 
'Trigeminy 1' type pathological heartbeats that were taken from a second 
patient. Fig. lie shows isolated 'PVC type pathological heartbeats that were 
taken from a third patient, and Fig. lid shows 'Bigeminy' type pathological 
heartbeats that were taken from a third patient. The problematic heartbeats 
were originally marked by red color, which could not be visible in the figures. 
Therefore, the red color has been replaced by horizontal lines to indicate the 
corresponding segments to which the description refers. 

Figs. He to llh show the different pathological heartbeats morphologies as 
detected by the system, per patient. For example, Fig. He shows a typical 
pathological heartbeat of 'patient (a)' (i.e. shown in Fig. 11a), which has been 
time-wise stretched in order to allow a close examination of its nature. 
Likewise, Figs, llf to llh refer to patient (b), (c) and (d), respectively. 

Fig. 12 schematically illustrates the process of detecting a large accidental peak 
that was superimposed on heartbeats, according to a preferred embodiment of 
the present invention. As described before, the peaks of the R-waves are 
detected by employing the WTA on the ECG signal. However, the latter 

detection is based on using a global maximum (i.e. per ECG segment) as 
reference, which could lead to heartbeat misdetection if there is an artifact with 
relatively large peak in the analyzed ECG segment. Therefore, such artifacts 
must be identi&ed and removed/discarded fi-om the ECG segment. Numeral 
reference 1201 is the acquiring part of the system. The sampled ECG signal 
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(1201a) is superimposed on a data picked-up by movement sensor 1201c, in 
order to establish a Base-Line (BL) reference. Alternatively, at the absence of a 
movement sensor, the EGG signal is lowpass-fQtered' (1201b), for obtaining the 
required BL reference signal (1208). In connection with BL signal (1202), a 
'Local Maximum Peak' process (1204) is utilized, for automatically detecting the 
most significant Local 'Maximum' peaks. Whenever a peak-point (not shown) is 
detected in the EGG signal (1201a), a corresponding movement (1209) is 
searched for in BL reference signal 1208. If such a movement is detected, the 
starting and ending instants of which being marked as timeg and time^ (in 1208), 
respectively, and the EGG peak 'falls' between said starting and ending 
instants, the corresponding EGG peak, together with its corresponding 
surroundings, is discarded and ignored (1205), because it is suspected as an 
artifact. The Local 'Maximum-peak' process utiKzes as a reference the Global 
Maximum of the entire new signal (1203), to which the local Maximum is 
compared. After discarding the artifact global maximum, the next global 
maximum is continued to be searched for within this segment, which is not 
suspected as artifact. Whenever such global maximum is found (i.e., that is not 
an artifact), it is utilized as a reference for the next EGG peaks. Nevertheless, 
the global maximums suspected as artifacts are saved, along with the rest of the 
local peaks, for later manipulations. 

If several consecutive local maximomis were found within an EGG segment, 
which were not suspected as artifacts (1206), the system refers to the largest 
maximum peak detected (i.e., within the EGG segment) as the global maximTun 
for this particular EGG segment (1207). 

Fig. 13 schematically illustrates adaptive detection of low frequency 
movements, according to a preferred embodiment of the present invention. The 
most problematic types of artifacts are those in the spectral range of the EGG 
signal, caused by patient movements, cable movements, etc. Artifacts having a 
frequency content that is similar to that of the actual EGG signal being 
recorded evidently can interfere with the correct interpretation of the EGG 



wo 03/057025 PCT/IL03/00005 
-56- 

signal. In order to detect this type of artifacts, a base line is extracted from the 
ECG signal, which represents the status of the monitoring hardware and the 
inherent physiological status of the monitored person, and an adaptive process 
is utilized for detecting changes in the base line. 

The system finds the 'normal' R-peaks (1301), and since their magnitude is 
likely to vary in time, the variance is averaged and the mean value is used as a 
reference (1302). The adaptive movement detection checks if there is a change 
in the variance (1305) of the base line in the environment of a heartbeat that is 
suspected as pathological (1303), according to the reference that was calculated 
adaptively according to the environment of the normal heartbeats (1302). If the 
ratio between the current variance to Ref. (1302), is larger than a threshold 
value (1305), it indicates that a movement has been detected. Otherwise, the 
system continues to find R-peaks (1301) and to update the Ref. (1302). 

Fig. 14 schematically illustrates adaptive detection of high frequency EMG and 
noise signals that are superimposed on viable heartbeats, according to a 
preferred embodiment of the present invention. This type of noise is 
characterized by having high frequency components caused by external noise 
and by muscle activity of the monitored person. In order to detect this type of 
noise, an adaptive detection process was developed, which utihzes the high 
scale of the Wavelet Transform Process (WTA) (1401). Using the WTA aUows 
analyzing the energy of the surroundings of individual R-waves, and 
determining a corresponding Ref. (1402). Than, the energy of the R-wave is 
analyzed and a corresponding reference (i.e., REnergy> 1403) is calculated. The 
ratio between the two references (i.e., Rg„e^lRef.) is calculated (1404), and if 
this ratio is larger than a threshold value, it indicates that the current 
heartbeat is noiseless (and the system continues by analyzing the next R, since 
in a noiseless heartbeat, as reflected in the ECG signal, the R-wave energy is 
much stronger than its surroundings. Otherwise, the current heartbeat is 
considered as noise (1405). 
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The analysis of the energies are analyzed, and the corresponding references 
(i.e., Ref. and Rg^rgy) recalculated, per R-wave and its surroundings, thereby 
allowing the system to respond very quickly to noises and EMG events, i.e., 
essentially immediately after their occurrences. The fast response, as described 
above, allows the system to utilize essentially every monitored heartbeat (i.e., 
BUTs) for enhancing the accuracy at which a person's cardiac condition is 
analyzed. 

Fig. 15 illustrates a general block diagram of the automatic Sleep Staging 
System, according to a preferred embodiment of the invention. The disclosed 
system is configured and programmed to identify the six known sleep stages 
(i.e. awakening, stages 1, 2, 3 and 4, and dreaming). Filter block, indicated by 
reference numeral 1501, comprises four Finite Impulse Response (FIR) filters 
(not shown), which eventually produce three types of signals. The first type is 
an EMG signal, which is typically superimposed on the EEG signal as 
interference. The second type is a Slow Frequency EEG (SFEEG), and the third 
signal is a Fast Frequency EEG (FFEEG). 

Adaptive segmentation block, indicated by reference numeral 1502, is an 
process that is based on an process known as the Generalized Likelihood Ratio 
(GLR) Test (a reference could be made to 
http://www.cnmat.berkeley.edu/~tristan/Report/node3.html), which is the basic 
statistical approach for sequential signal segmentation. The mathematical basis 
of the GLR will not be described here since it is known to those skilled in the 
art. The task of the GLR-based process is to partition the resulting filtered 
signals (i.e. SFEEG, FFEEG and EMG) iato quasi-stationary epochs, for 
allowing identifying and extracting features (1503) that are contained therein. 
These features are: (1) Mean value, (2) Variance, (3) Skewness, (4) Duration, (5) 
Relative Energy in Delta Frequency Band (1-3.5 Hz), (6) Relative Energy in 
Theta Frequency Band (4-7 Hz), (7) Relative Energy in Alpha Frequency Band 
(7.5-12 Hz), (8) Relative Energy in Sigma Frequency Band (12.5-15 Hz); and (9) 
Relative Energy in Beta Frequency Band (15-35 Hz). Features Nos. (1) to (4) 
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are time-wise features, and features Nos. (5) to (9) are spectrum-wise featxires. 
All of the nine features are extracted from the (filtered) EEG signal, but only 
features Nos. (1) to (4) are affected by the EMG signal. 

Identification of the above-specified features allows detecting, by utilizing the 
block indicated by reference numeral 1504, the following events: (1) K- 
complexes, (2) Body movements, (3) Electrode movements, (4) Rapid eye 
movements, (5) Eye blinks, (6) Spindles; and (7) Alpha intrusions. Events Nos. 
(1) to (5) are time transients and are easily extracted from the ongoing EEG 
signal due to their distinguishable energy, skewness and duration in relation to 
said EEG signal. The time transient events (i.e. events 1 to 5) are classified 
using their new environment frequency behavior and their time domain 
morphology. Events Nos. (6) and (7) are frequency transients and are detected 
by measuring their duration and by analyzing their spectral behavior. 

In order to obtain only pure features of ongoing sleep stages from the EEG 
signal, an Adaptive Frequency Activity Detection (AFAD) modide is utihzed, 
indicated by reference nimieral 1505, which is based on the Hidden Markov 
Model (HMM), which is also known as the "Forward-Backward Process". HMM 
models are utilized for extracting signal patterns of interest, by characterizing 
their frequency variations. In the present invention, these models are utihzed 
for optimizing a patient-learning process hy analyzing three major groups of 
fi:equency activities (i.e., High Frequency, Mixed Frequency and Low 
Frequency). The functionahty of the HMM model will not be discussed since it 
is known to those skilled in the art. A notion, regarding HMM models, coidd be 
fotmd in the article "Hidden Markov Model (EMM) Toolbox", which is given in 
the Internet web site http://www.cs.berkeley.edu/~murphyk/Bayes/hmm.html. 

Another model used for learning individual patients is the Abstract State 
Machine (AMN) model, a notion of which could be found in a book called 
"Software Engineering with B", by J.B. Wordsworth, Published August 1996 by 



wo 03/057025 PCT/IL03/00005 
-59- 

Addison Wesley Higher Education. For additional details visit the Internet web 
site http://www. pear8onptg.com/book_detaiy 0,3771,0201403560,00.html). 

After detecting the time transient events, their features are subtracted from the 
features signals, resulting in pure features sleep EEG only. The spectral 
features 5, 6, 7 and 9 (relative Energy in Delta Frequency Band, Theta 
Frequency Band, Alpha Frequency Band and in Beta Frequency Band, 
respectively) are iaput into two HMM model types, as illustrated in Fig. 18, for 
detecting frequency activities. 

Fig. 16 schematically illustrates extraction of the SFEEG, FFEEG and the EMG 
signals from the raw EEG signal, according to a preferred embodiment of the 
invention. Filter 1601, being FIR 64-order Band Pass Filter (BPF) 1-35 Hz, 
essentially removes the EMG signal ffom the raw EEG signal and cancels drifts 
in the base line of the EEG signal. The filtered EEG signal is, then, decomposed 
by two filters - filter 1603, being FIR 64-order BPF 1-7 Hz, and filter 1604, 
being FIR 64-order BPF 7-35 Hz, into SFEEG signal and FFEEG signal, 
respectively. PQter 1602, being FIR 64-order BPF 30-35 Hz, extracts the EMG 
signal for further analysis. 

Fig. 17 is a block diagram illustrating adaptive identification of the boundaries 
of duration-variable segments in the SFEEG signal, in the FFEEG signal and 
in the EMG signal, according to a preferred embodiment of the invention. 
Adaptive segmentation process 1701 is employed on the SFEEG signal, for 
identifying segments that are associated with corresponding sleepiag stages. 
Likewise, adaptive segmentation processes 1702 and 1703 are employed on the 
FFEEG signal and on the EMG signal, respectively. The SFEEG and FFEEG 
boundaries are then combined (1704). After identifying the above-mentioned 
boundaries, the corresponding featvires and special events are extracted (Fig. 
15, reference numerals 1503 and 1504, respectively). 
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Fig. 18 schematically illustrates detection of frequency activities by utilizing a 
three-state HMM model (a) and a two-state HMM model (b), according to a 
preferred embodiment of the present invention. The incentive for utilizing the 
HMM models is the ongoing sleep EEG signal being a Markovian process, 
namely, the current status of the EEG signal is largely affected by its past 
behavior. Therefore, only several particular transfers, from one stage to 
another, are allowed from statistical point of view. According to the invention, 
the frequency spectrum involved in sleep stages is divided into three groups of 
frequency activity, each group indicating a different sleep stage: (1) Fast 
Frequency Activity (FFA), (2) Mixed Frequency Activity (MFA); and (3) Slow 
Frequency Activity (SFA). The latter groups being indicated by reference 
numerals '1', '2', and '3' (see Fig. 18a). The goal in utilizing the HMM models 
is to decide which model (i.e. (a) or (b)) better describes the sleep status at a 
given time. Model (a) indicates that the patient was sleeping deeply, and model 
(b) indicates that the patient was having a 'Ught' sleep. Each stage is 
characterized by a 4-diniensional probabiUty Gaussian. Passing from one type 
of activity (i.e. sleep stage) to another is characterized by a corresponding 
probability. For example, the probability of passing from state 1 (i.e., FFA) to 
state 2 (i.e., MFA) is higher than the probability of passing from state 1 to state 
3 (i.e., SFA). 

The decision criteria, according to which the three-state model, or the two-state 
model, is chosen, is the relative energy (in percentage) of the delta-wave. If the 
latter energy is larger than 40%, the three-state model is employed. Otherwise, 
the two-state model is employed. Prior to deciding which model is valid, the 
system is first initialized by assigning it a set of predetermined global initial 
conditions (that are input into these models), as well as the patient's feature 
matrix. The next step is employing the BAUM WELCH training 

process, the task of which is to allow the system to converge from the initial 
statistical properties (i.e., being provided by a clustering process based on a set 
of fuzzy logic rules) of each model, into the patient's self-parameters. The 
BAUM WELCH training process outputs a matrix describing the probability of 
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a patient being ia each sleep stage essentially at any given time. The 
mathematical background of this process will not be discussed, since it is known 
to those skilled ia the art. For example, a reference could be made to the article 
"Derivation of the Baum-Welch process for HMMs", which could be foiind in the 
Internet web site 

http://www.speech.sri.com/people/anand/771/html/node26.html. 

The resulting probability matrix is averaged in time by frames of 30 seconds 
while avoiding overlapping, and entered, together with the time-wise events, 
into the R&K operator, in order to obtain a final decision regarding each sleep 
stage (Fig. 15: 1506). The R&K operator will not be discussed here, as it is 
known to those skilled in the art. A notion, regarding R&K rules, could be found 
in the article "Objectives of the SIESTA project", which could be found in the 
Internet web site http://www.ai.univie.ac.at/oefai/nn/siesta/objectives.html. 

Fig. 19 (prior art) shows an example of two consecutive quasi-stationary 
segments that were derived from typical EEG signal. In order to decide whether 
quasi-stationary segments X and Y belong to the same stationary process, the 
resemblance between' their respective statistical properties has to be 
determined. The greater the resemblance, the more likely that segments X and 
Y belong to the same stationary process. The extent of the resemblance 
X(0<x<l)is determined by utihzing the GLR process. 

Fig. 20 schematically illustrates the principle of sequential segmentation 
process for identifying quasi-segments in EEG signal, according to a preferred 
embodiment of the present invention. First, a search process is employed, for 
identifying a new current potential (temporal) boundary (2001) between two 

consecutive segments, the latter process being iterative process. Whenever a 
new boundary is found, a second process is employed (2003), for optimizing the 
position/location of the potential boundary. After the current boundary's 
location is optimized (2004), the next current boundary is searched for. Each 
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two consecutive detected boundaries define the duration of the corresponding 
segment (i.e., the segment 'locked' between these boundaries). 

Fig. 21 schematically illustrates the principle of identifying quasi-segments in 
BEG signal, according to a preferred embodiment of the present invention. The 
decision process that is utilized for detecting a new boundary at an arbitrary 
point 'n' (Sb<n<Se+L-l) in a sampled EEG signal is performed by utilizing a 
Sliding Test Window (STW), marked 'S', and a Growing Reference Window 
(GRW), marked 'R'. Sequences {R} and {S}, each of which includes 
corresponding samples of the raw EEG signal, are compared to one another by 
employing the GLR test. If the (statistical) distance between the latter 
sequences exceeds a Predetermined Threshold Value (PTV), the GRW and the 
STW are gradually shifted to the right while maintaining their initial width. 
The right-shift continues until the GLR distance is smaller than the PTV. At 
this stage, the optimal boundary (i.e., between two consecutive distinguishable 
segments) location is assiuned to be between Sb and Se, and a process for 
determining the optimized (i.e., more accurate) boundary location is employed. 
However, if the GLR distance initially does not exceed the PTV, Se is 
incremented causing the GRW to grow and the STW to slide to the right, until 
the GLR distance exceeds the PTV. At this stage, the optimal boundary location 
is assumed to be between Se and Se+L-1, and the process for determining the 
optimized boundary location is employed. The above-described process, for 
determining a new boundary, is repeated until all the boundaries are identified. 

Figs. 22a and 22b are exemplary results showing the envelope of an effort 
signal and the corresponding segmentation of the effort signal, according to a 
preferred embodiment of the present invention. In Fig. 22a, the firequency and 

amplitude of effort signal 2201 change according to the Respiratory rhythm and 
the physiological status of the monitored person (not shown). Each peak-point 
in effort signal 2201 is identified, and its amplitude is represented by a 
corresponding point. For example, the peak-points 2202a, 2202b and 2202c are 
represented by Amphtude Points (APs) 2204a, 2204b and 2204c, respectively. 
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Envelope 2203 comprises all the APs. The next step is smoothing the signal 
according to fuzzy logic rules and finding APs that are global maxima, after 
which their corresponding instants wiU be utilized as time-boundaries for 
segmenting effort signal 2203. 

In Fig. 22b, APs 2205 to 2209 are global maxima and are indicated by 
corresponding vertical lines, such as line 2211. AP 2210 is also a global 
maximum. However, since it is too close (according to a certain criteria) to AP 
2207 and smaller than AP 2207, it is not utDized as time-boundary. 
Accordingly, effort signal 2201 is segmented into segments SI to S6, and each 
segment is evaluated. 

Fig. 23 shows one exemplary segment of the effort signal, according to a 
preferred embodiment of the present invention. Segment 2300 is evaluated in 
order to decide whether it reflects an apnea or hypopnea. Accordingly, the 
magnitude of the mean value of points close to minimal point 2302 is compared 
to the magnitude of the mean value of points close to the magnitude of maximal 
point 2301, for evaluating the probability of a magnitude decrease of more than 
20%. Another criteria is implemented, according to which it is decided whether 
the trend of the segment is Up, Down or Even (i.e., 'flat')- If the segment is flat, 
a corresponding PWW (not shown) is utilized by the FLA, for evaluating the 
correctness of the assumption e.g. if a respiratory event (hypopnea or apnea) 
occurred. 

Figs. 24a and 24b (prior art) show exemplary correlation between an effort 
signal that contains segments in which magnitude decrease of more than 20% 
and for duration of more than 10 seconds were detected, and Oxygen 
Saturation. As can be seen in Fig. 24a, significant Respiratory (i.e., effort) 
deterioration in, e.g., sections A and B (i.e., apneas) caused significant 
deterioration in the Oxygen Saturation, as reflected in sections C and D, 
respectively. 
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Fig. 25 schematically illustrates adjoining Downwards-trend segments and 
Upwards-trend segments, according to a preferred embodiment of the present 
invention. After finding parabola 2501, which approximates the envelope of 
segment 2500a, and parabola 2502, which approximates the envelope of 
segment 2500b, the Segment Combining process (SCA) is employed, after which 
a new segment 2500c is obtained. The 20% reduction feature is searched for in 
the new (grand) segment 2500c. 

Fig. 26 schematically illustrates exemplary combinations of different events, on 
which the Decision-Making Process (DMA) is employed, according to a preferred 
embodiment of the present invention. The chest signal reflects the effort signal. 
The ' mark indicates the beginning of an event, and the '© ' mark indicates 
its ending. In this example, four apnea events have been identified, i.e., lA, IB, 
IC and ID. The duration (i.e., time between '• ' and corresponding ' © ') of each 
of the apnea events lA to ID is random. Since a decrease in the abdominal 
eJEfort is expected to overlap the chest effort, a PWW that is associated with the 
abdominal event is located in relation to the apnea event in a way that the '1.0 
probabiHty' (horizontal) line of the abdominal event starts and ends essentially 
at the starting point (i.e., •) and ending point (i.e., ®) of the apnea, 
respectively. For example, the starting and ending points of abdominal event 
2D (i.e., 2601 and 2602) essentially overlap the starting and ending points of 
chest event ID, respectively. Accordingly, abdominal event 2D is assigned the 
largest weight (i.e., probability value 1). Abdominal event 2A is assigned a 
smaller weight, because its starting point 2603 is almost outside the limits of 
the PWW window 2605. Likewise, hmb events 3A and SB, Oximeter events 4A 
and 4B will be assigned the maximum weii^ht. 

Eegarding chest event IG, an Oximeter event 2606 has been detected. However, 
since it is completely outside PWW 2607, it is assigned a zero weight (value). 
Namely, Oximeter event 2606 can not be logically related to any other event. 
Consequently, Oximeter event 2606 is ignored and therefore all D event. 
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Fig. 27 (prior art) is a graph illustrating resulting apnea events according to 
human scorer. Reference numerals 2701a to 2701d represent aggregations of 
apnea events. Because conventional systems are incapable of deciding the 
probability of occurrence of an apnea event, each event is assigned either a 
unity value or a zero value. 

Fig. 28 is a graph illustrating resulting apnea events, as affected by other 
factors, according to a preferred embodiment of the present invention. 
Reference numeral 2801a is an aggregation of several events that are suspected 
as apneas. In these aggregated events, additional two tj^jes of events were also 
considered in order to evaluate the probability of an event being an apnea. One 
type of these events is the Oxygen saturation level in the blood (i.e., as 
measured by the Oximeter channel), and the second type is the probability that 
the monitored person was sleeping (i.e., measured by the EEG signal). For 
example, at time tl the resulting probability for an apnea is 0.3. A low- 
probabihty apnea may be obtained as a result of having: (1) very low probabihty 
value that the detected corresponding decrease in the Oxygen saturation level is 
associated with an apnea, and/or (2) very low probability value that the 
monitored person was ftdly asleep. At time t2, however, the resulting 
probability of an apnea is much higher (i.e., 0.82), because at that time the 
probabihty that the monitored person was sleeping was relatively high, and/or 
the corresponding decrease in Oxygen saturation level was more closely related 
to the event that was suspected as apnea. 

Reference numeral 2801b refers to another aggregation of events. In these 
aggregated events, only one additional type of event was also considered in 
order to evaluate the probability of an event being an apnea, namely, the 

Oxygen saturation level in the blood. In general, in this aggregation, apnea 
events were assigned relatively high probabihty values. 
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In general, the more event types (e.g., Oxygen saturation. Sleep status) are 
considered when evaluating the probability of an apnea, the more 'realistic' is 
the decision regarding such apnea. 

Fig. 29 schematically illustrates the functionality and layout of the system, 
according to a preferred embodiment of the present invention. Channels 1 to 7 
are associated with ECG signal, EEG signal, EMG signal and Respiratory 
signals - abdominal effort signal, chest effort signal and Oxygen saturation 
signal, respectively. Samples of these signals (not shown) are input into a 
computer, which contains the software processes 2901 that are required for 
taking corresponding decisions based on the processed samples. The system is 
capable of learning individual heartbeat behavior, as reflected by the 
corresponding ECG signal, as well as sleep behavior, as reflected in the 
corresponding EEG signal. Moreover, the system is also capable of learning 
movements that characterize the monitored person, thereby allowing 
distinguishing normal movements jfrom random movements, which are usually 
regarded as noise. The learning capabilities of the disclosed system allow 
relating to each other the ECG, EEG, EMG and the Respiratory signals (2901a 
to 2901c), which allow to obtain more accurate and enhanced decisions. For 
example, regarding Respiratory disorders, an accurate decision may be obtained 
after considering also the corresponding sleep stage (i.e., EEG signal) of the 
person. 

Storage array 2903 is utihzed in general for storing several types of data: (1) 
samples of the raw signals (e.g., ECG), (2) initial parameters, for allowing the 
system to start its operation whenever a person is connected to the system, (3) 
processed data, which is associated with the system learning process of a 
monitored person. 

Displays Dl to Dn may display various specific physiological data of interest, on 
request, according to application layer 2902, such as time-wise events and/or 
frequency-wise events. For example, display Dl may introduce sleep stages, D2 
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heartbeats, D3 Eespiratory signal, etc. If a clinician wishes so, a display may 
introduce him with raw data of desired type(s) (e.g., ECG). 

Software 2901, Dl to Dn and application layer 2902 may reside within same 
computer or in different locations, in which case they may communicate with 
each other over a data network 2904, such as the Internet, LAN etc. 

Fig. 30 schematically illustrates a CeUular-phone-based Remote Monitoring 
System (CRMS), according to a preferred embodiment of the present invention. 
Stripe 3001, which includes, e.g., electrodes 3001a and 3001b, is placed on a 
monitored person (not shown). Electrodes 3001a and 3001b include transmitters 
(not shown) for transmitting, over wireless communication channel 3002a, 
relevant signals (e.g., ECG, EEG) to ceUiilar phone 3002, which is located close 
enough to the transmitting antenna(s), to allow good reception of these signals. 
Alternatively, electrodes 3001a and 3001b are connected to interface 3003, 
which transmits these signals to cellular phone 3002 on wireless 
communication channel 3002b. Transmissions over wireless communication 
channels 3002a and 3002b may be carried out by utilizing the Wireless 
AppHcation Protocol (WAP) standard. 

Cellular phone 3002 contains software, the task of which is to carry out a 
prehminary classification process of the signals, i.e., for identify typical and 
irregular shapes in these signals, and transmit the latter signals to cellular 
base-station 3004, in which a dedicated server 3005 resides. Server 3005 
contains software, the task of which is to carry out a fuU evaluation of the 
received signals, essentially as described above. Whenever a decision is made in 
server 3005, regarding heartbeat (or heartbeat rate) disorder, and problematic 
segments (i.e., of the types described above), the corresponding data is 
transmitted from base-station 3004 to clinician's cellular phone 3006, and 
displayed on its display in the form of a graph or otherwise. 
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There are several drugs that only a qualified physician can authorize their 
usage. However, a qualified physician is rarely present at a scene wherein a 
person urgently needs such a drug, and a considerable amount of time may be 
lost until the person is brought to such qualified physician (e.g., in a hospital). 
In some cases, the time lost may be critical to the life of the person. 

The CRMS system solves the latter described problem by allowing paramedical 
staff to administer drugs to a person while taking orders fi:om a physician that 
is located far away firom the paramedical staff. According to the invention, upon 
arrival at a scene, wherein a person is found to be in medical critical situation, 
an, e.g., ambulance staff attaches the electrodes to the person, and the relevant 
data is received by a physician's cellular phone, essentially as described above. 
Upon receiving the latter data, the physician orders the paramedical staff to 
administer to the person a certain dosage of corresponding drug, while the 
person is being monitored. As the person's physiological status changes as a 
reaction to the drug being administered to him, and the data appearing on the 
physician cellular phone's display changes accordingly and 'on-line', the 
physician may give the paramedical staff correcting orders (i.e., by utOizing a 
different, or same, cellular phone), until the status of the person stabilizes or 
improves, whichever the physician's decision is. 

Of course, the cellular phone system coidd be replace by other types of (data) 
networks, such as the Internet, a Local Area Network (LAN), a Wide Area 
Network (WAN) etc. 

Fig. 31 schematically illustrates utilizing the Automatic Sleep Stage (ASS) and 
Respiratory system with combination with a Continuous Positive Airway 

Pressure (CPAP) system, according to a preferred embodiment of the present 
invention. The CPAP system is the most common treatment for sleep apnea. 
The patient wears mask 3101 over his nose during his sleep, for allowing 
controllable CPAP system 3104 to force air through his nasal cavity to provide a 
pneumatic spHnt that prevents the upper airway fi:om obstructing. The air 
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pressure is adjusted so that it is just enough to prevent the throat from 
collapsing during sleep. Electrode 3102 is mounted in mask 3101 in such a way, 
that it is tightly attached to the relevant spot on the face of the monitored 
person to allow picking up the required signal(s). Electrode 3102 transmits 
(3105) signal(s) to ASS & Respiratory system 3103, which processes the 
signal(s) essentially at the way described above. Whenever required, system 

3103 sends a signal to CPAP system, for causing it to change the pressure of the 
air that is provided to mask 3101. The latter signal is associated with the 
cxurrent sleep stage and Respiratory condition as detected by system 3103. 

Conventional CPAP systems (i.e., that do not include ASS systems), have 
several drawbacks: (1) apnea events may return whenever the CPAP is used 
improperly, which is normally the case because conventional CPAP systems are 
poorly controlled, or not controlled at all, by physiological feedback. Therefore, 
these systems output a constant and continuous pressure, or vary the pressure 
to coincide with the person's breathing pattern, or start with low pressure, after 
which the pressure is gradually increased. Pressure is constantly provided to 
the person's nasal cavity causii^ to side effects, such as nasal and facial skin 
irritations and drying, abdominal bloating, mask leaks, sore eyes, breathing 
interference and headaches. From technical point of view, conventional CPAP 
systems are inefficient, because they supply, at different periods, air that is not 
beneficial to the patient. 

Utilizing ASS & Respiratory system 3103 ia combination with CPAP system 

3104 is advantageous because system 8103 is capable to foresee an apnea event, 
in which case it causes the CPAP to provide the patient with the exact air 
pressure at the right time. In cases where system 3103 does not foresee an 

apnea, no pressure is employed by the CPAP system. Therefore, affihation of 
system 3103 with CPAP systems greatly improves the performance and 
efficiency of the CPAP systems, and the side effects are significantly reduced. 
Due to system 3103 adaptation capabihty, to adjust (i.e., by learning) itself to 
the monitored patient, the CPAP also adapts itself to the monitored patient. 
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Monitoring xinit 3107 allows local, or remote, supervision of a clinician, which 
might be passive, due to the 'automatic nature' of the system, or active, due to, 
e.g., emergency state. 

Figs. 32a and 32b illustrates exemplary results of the adaptive segmentation of 
EEG and EOG signals, respectively, according to a preferred embodiment of the 
present invention. In Fig. 32a, vertical liaes 3201a to 3202i indicate the 
boundaries of the corresponding quasi-stationary segments, which were 
calculated by the adaptive process as described above. As might be appreciated, 
each quasi-stationary segment has a unique shape, which distinguishes it from 
its adjacent quasi-stationary segments. For example, segment SI differs in 
shape from segments S2 and S3. Likewise in Fig. 32b, wherein EOG quasi- 
stationary segments were foxmd in the same way as the EEG quasi-stationary 
segments. 

Fig. 33 illustrates an example for associating quasi-stationary segments with 
one of the three frequency groups, according to a preferred embodiment of the 
present invention. After extracting the time and frequency features from the 
quasi-stationary segments, each quasi-stationary segment is associated with a 
corresponding frequency group. For example, quasi-stationary segments 
indicated by 'H' were found to contain relatively high-frequency content. 
Consequently, these segments wotdd be associated with the FFG. Likewise, 
quasi-stationary segments indicated by 'M' would be associated with the MFG, 
and quasi-stationary segments indicated by 'L' would be associated with the 
SFG. 

Fig. 34 illustrates detection of special events, according to a preferred 
embodiment of the present invention. Special events are identified by their low, 
or high, frequency content and by their temporal shape. Segments marked with 
'L' are low-frequency events, and segments marked with 'H' are high-frequency 
events. Special events are not utihzed in the cluster analysis. 
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Figs. 35a to 35cl illustrate exemplary classification of special events, according 
to a preferred embodiment of the present invention. In Fig. 35a, typical K- 
complex is shown (marked with 'K'). In Fig. 35b, typical spindle is shown 
(marked with 'S'). In Fig. 35c, typical eye-bMnk is shown (marked with 'EB')- In 
Fig. 35d, tjTpical body movement is shown (marked with 'BM'). 

Figs. 36 illustrate exemplary REM events, according to a preferred erabodiment 
of the present invention. REM events are detected only if the ASS system is 
configured to handle EOG signal, typical REM events are marked (in Fig. 36) 
with 'R'. 

Fig. 37 illustrates exemplary clusters, according to a preferred embodiment of 
the present invention. Graph 3700 shows typical clustering of quasi-stationary 
segments features as were detected in the whole night. Axis 3701 and axis 3702 
describe the relative energy (between 0% and 100%) of Delta-wave or Alpha- 
wave, respectively. Axis 3703 describes the total energy of the corresponding 
quasi-stationary segment. White colored cluster 3704 describes every quasi- 
stationary segment that was associated with the FFG group. Gray colored 
cluster 3705 describes every quasi-stationary segment that was associated with 
the MFG group, and the black colored cluster 3706 describes every quasi- 
stationary segment that was associated with the SFG group. 

Each cluster (i.e., 3704, 3705 and 3706) is characterized by a Gaussian 
Probability Function (GPF, not shown). Each fi-equency group (i.e., cluster) is 
averaged per 30-second time-interval, resulting in a membership value (i.e., 
between 0.0 and 1.0). Accordingly, each 30-second time-interval is represented 
by corresponding three membership values, each of which represents the 
probabiKty of a person being in a corresponding sleep stage. The actual sleep 
stage (i.e., per specific time -interval) is obtained by employing the R&K set of 
rules on the corresponding membership values and on the relevant detected 
special events. Fig. 38 depicts the resulting hypnogram, which describes the 
relevant sleep stages. 
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Figs. 39a and 39b illustrate exemplary detection of low energy portions in the 
EMG signal, according to a preferred embodiment of the present invention. 
Sections denoted LE represent the low energy portions of the EMG signal. In 
Fig. 39b, reference numerals 3901 to 3907 represent the membership values of 
the EMG signal. 

Figs. 40a to 40c illustrate typical Awakening, Arousal and Micro-Arousal 
stages, according to a preferred embodiment of the present invention., 'A' (In 
Fig. 40a), 'B' (bx Fig. 40b) and 'C (In Fig. 40c) denote the Awakening, Arousal 
and Micro- Arousal stages, respectively. 

Fig. 41 is a 2-dimensional representation of relative activity of (the three) EEG 
frequency groups, according to a preferred embodiment of the present invention. 
Fig. 41 is, in fact, a 2-dimensional (2D) complimentary representation of the 3- 
dimensional representation of the exemplary EEG frequency groups' activity 
shown in Fig, 87. 

The original graphs, the pattern of which are shown in Fig. 41, have different 
colors for allowing distinguishing between different EEG frequency groups. 
However, since Fig. 41 contains only a black and white version of the latter 
pattern, the actual colors have been replaced by corresponding references (i.e., 
'BL', 'GR' and 'RD'). Accordingly, 'BL', 'GR' and 'RD' denote the relative 
activity percentage of high frequency EEG, mixed frequency EEG and low 
frequency activity, respectively, as a function of time. 

In Fig. 41, the "X" axis denotes the time axis, and the "Y" axis denotes the 
relative average activity (in percentage) of a frequency group (with respect to 
other frequency groups). Each point of each one of the three exemplary graphs 
was obtained by averaging the activity level of the respective frequency group 
over a corresponding 30-seconds time epoch. However, time epochs of lengths 
other than 30 seconds are optional. For example, points 'A' and 'B, which are 
related to graphs RD and GR, respectively, were obtained by averaging the 
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activity of the respective frequency groups over time epoch 7^ (preferably 
having duration of 30 seconds). Accordingly, points A and B were found to have 
87% and 13% of the overall activity related to time epoch . The remaining 
frequency group (i.e., represented by the 'BL' graph) showed no activity at all 
dxiring time epoch Tg and, therefore, point 'C has 0% relative activity. 

The advantage of the 2D representation lays in the fact that the 2D 
representation shows, simultaneously, several types of brain activities that 
could be represented by frequency groups, while showing the relative activity of 
each one of the related brain activity tj^es. Furthermore, the 2D representation 
is advantageous from the clinical point of view, because the 2D representation 
offers an efficient way for identifying a patient's specific sleep pattern, 
including the dynamics of his sleep, which is important in diagnosis of sleep 
disorders, such as narcolepsy, insomnia and alpha intrusion. 

Fig. 42 schematically illustrates a 'one-Lab-Server' system according to one 
most preferred embodiment of the present invention. The system comprises 'm' 
Gateways (i.e., gateways 4201/lA to 4201/mA), each of which could be, e.g., a 
personal computer (PC) residing within a patient's residence or within a private 
or public cHnic, 'n' User Workstations (i.e., user workstations lA to nA), and 
Lab Server 4202A. To each one of (Gateways 4201/lA to 4201/mA could be 
connected at least one monitor, such as monitor 4204, and/or at least one 
treatment device, such as treatment device 4205. A monitor is an apparatus 
capable of acquiring biomedical signals from a person that is connected thereto 
by utilizing corresponding electrodes, or sensors, and a treatment device is an 
apparatus capable of accepting external control signals (e.g., from local or 
remote computer) and administering drugs into a person connected thereto, 
according to the control signals. The Gateway is responsible for converting the 
acquired biomedical signals into corresponduig formatted raw digital data, 
recording and forwarding the formatted raw digital data to central Lab Server 
A (4202A), via a data network (not shown), where the formatted raw digital 
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data is processed, analyzed and classified, according to the method disclosed in 
the current invention. Each one of Users Workstations lA to nA allows a user, 
e.g., a physician, an access to Lab Server 4202A, thereby allowing the 
user/physician retrieving and reviewing selected raw and/or analyzed data 
(such as EEG, respiratory and EGG signals) that is related to certain tinie(s) 
and/or aspect(s) of the physiological condition of a patient that is, or was, 
connected to any of the Gateways. For example, a physician may retrieve data 
from Lab Server 4202A for reviewing the hypnogram, sleeping statistics, and 
special/meaningful events, which are related to heart and respiratory condition 
of a specific patient. The interaction between a physician and central Lab 
Server 4202a may be in the form of inquiries and/or requests. For example, a 
physician may request the 'Patient-Specific Universal HeartBeat' (or other 
physiological/event-related data) of specific patient, regardless of the Gateway 
to which the patient is, or was, connected via a corresponding monitor. 
According to another example, the physician may request data that is related to 
the heart condition of all patients that are, or were, treated by the same 
physician, for, e.g., statistics purpose. According to still another example, the 
physician may request data that is related to every physiological behavior, 
regardless who the patient is, that matches a same specified pattern, for, e.g., 
research purpose. Every request is forwarded from the physician workstation 
(e.g., Workstation No. 2A) to Lab Server 4202A. Lab Server 4202A responds to 
the request by searching for the requested data in a database that could reside, 
e.g., within Lab Server 4202A. After being found by Lab Server 4202A, Lab 
Server 4202A generates an XML / HTML file that contains the requested data, 
and forwards the XML / HTML file to the requesting Workstation (e.g., to 
Workstation 2A). The limit on the actual number of Gateways, which could be 
handled by a Lab Server, depends on the capacity of the Lab Server. Currently, 
one Lab Server is capable of handling about 40 Gateways. 

A bi-directional communication channel, between Lab server 4202A and each 
one of the Gateways lA to mA, over a data network (not shown), allows 
biomedical data to be forwarded from Gateways lA to mA to Lab server 4202A, 
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and control signals to be forwarded from Lab Server 4202A to the respective 
Gateway, in order for the control signals to be relayed to monitors (e.g., for 
initialization of the monitors, setting their gain etc.) and treatment apparatuses 
(for, e.g., determining the rate flow of a drug and the timing for administering 
the drug into a patient). The control signals could be initiated by Lab server 
4202A and/or by a user, based on the displayed data in a user workstation. 
Fig. 43 schematically illustrates a system according to Einother preferred 
embodiment of the present invention. Basically, the system shown in Fig. 43 
may contain as many as required systems as shown in Fig. 42, in modular 
manner, and also a central server for maintaining coordination between aU 
components of the new system of Fig. 43. The system may comprise additional 
LAB Servers, such as Lab Servers 4202B and 4202C, each of which functions 
essentially in the same manner as described in connection with Lab Server 
4202A of Fig. 42. Central server 4301 functions to intermediate between all 
user workstations (i.e., user workstations A, B, C etc.) and all Lab Servers (i.e., 
Lab Servers A, B, C, etc.), and has URL and Internet network connection. 
Central server 4301 has also a central database 4303 that contains data related 
to each one of the databases residing within, or related to, every Lab Server, 
and in particular, database 4303 contains data related to users, their login 
information and authorization level information; gateways parameters; and 
patients' personal and cUnical data. Central Server 4301 functions to identify 
users, and to route or relay users to relevant Lab server(s), from which 
requested biomedical data is retrieved by the users. Central server 4301 also 
functions as the general manager of the multi-patient/multi-physician (i.e., 
user) system, thereby allowing automatic routing (connection) of users to the 
Lab Servers by the Central Server. Alternatively, the identification of users by 
the Central Server can be performed by requesting and receiving information 
from Authentication and Encryption Server 4302. 

Gateways are connected directly to the corresponding Lab Servers. For 
example, Gateway 3A is coimected to Lab Server 4202A, Gateway 7B is 
coimected to Lab Server 4202B, and Gateway 5C is connected to Lab Server 
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4202C. However, Gateways 3A, 7B and 5C are exemplary Gateways shown only 
for illustration purpose. In fact, to each one of the additional Lab Servers 
(4202B, 4202C, etc.) could be connected, via a data network, a pluxaUty of 
Gateways and user workstations, essentially in the same manner as described 
in connection with Lab Server 4202A of Fig. 42. 

Authentication server 4302 provides information for authentication of users and 
encryption keys for secured transmission. Whenever a user with certain 
permission level logs-in (via the respective user workstation), Central server 
4301 forwards the 'login data' to Authentication and Encryption server 4302, 
for authentication, which checks the permission level of the user, and allows the 
user to interact with the system (i.e., retrieving biomedical data and initiating 
control signals) according to this permission level. 

According to one aspect of the present invention, each one of the user 
workstations may forward a request for specific biomedical data to Central 
server 4801, after which Central server 4301 responds by searching in database 
4303 for a data that wUl be helpftd in routing the user to the corresponding Lab 
Server (i.e., the Lab Server containing the requested biomedical data), from 
which the user may retrieve the required medical data. 

Of course, individual one-Lab-Server systems may be configured to be self- 
contained. 

The above examples and description have of course been provided only for the 
purpose of illustration, and are not intended to limit the invention in any way. 
As will be appreciated by the skilled person, the invention can be carried out in 
a great variety of ways, emplosdng more than one technique from those 
described above, all without exceeding the scope of the invention. 
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1. Method for monitoring normal/abnormal physiological events in one or more 
patients by the analysis of biomedical signals of each patient, comprising: 

for each patient, 

a) aggregating one or more raw data streams, each of which representing a 
different type of biomedical signals of said patient, in a data storage; 

b) for each type of biomedical signal, performing adaptive segmentation of 
said raw data streams associated with said patient, and assigning, to each 
segment, individual attributes being represented by attribute values; 

c) determining an attribute domain, in which each segment being 
represented by a point that corresponds to the attribute values of said segment; 

d) for each tjrpe of biomedical signal or for any combination thereof, 
generating a set of clusters in said attribute domain, each of which consisting of 
a combination of points determined by their relative location to other points, by 
assigning a set of property values to each point, each property value 
corresponding to the degree of association of said point with one of the clusters; 

e) obtaining the probability of the order of appearance of each point, in time, 
according to its property value and updating the property value of each point in 
each cluster using said probability; 

f) repeating steps d) to e) above while in each time, varying the combination 
of points included in each cluster according to their most updated property value 
and by including points derived from adaptive segmentation of further 
aggregated raw data streams of said patient, thereby updating each cluster; 

g) associating at least one updated cluster with a normal/abnormal 
physiological state of said patient by using fiDrmer knowledge, regarding 
normal/abnormal physiological states of said patient and/or of a reference group 
of patients, that is represented as reference clusters in said domain; and 

h) individually characterizing said patient by identifying normal/abnormal 
physiological states, associated with one or more updated clusters, and utilizing 
the order of appearance of said normal/abnormal physiological states for 
predicting normal, and/or abnormal, physiological states of said patient. 
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2. A method according to claim 1, further comprising obtaining indications 
related to abnormal physiological events in the patient by displaying changes in 
normal physiological states of said patient as a function of time. 

3. A method according to claim 1, further comprising displajdng different stages 
of normal physiological states of the patient as a function of time. 

4. A method according to claim 1, further comprising storing at least a portion 
of the aggregated raw data streams in an accessible data storage. 

5. A method according to claim 1, further comprising retrieving the data 
segment that corresponds to an identified physiological disorder from the 
accessible data storage. 

6. A method according to claim 1, wherein the types of biomedical signals are 
selected from the group: 

- EEG signals; 

- EOG signals; 

- EKG signals; 

- Respiratory signals. 

7. A method according to claim 1, further comprising: 

a) Generating an individual data being a group of clusters and/or 
indications of normal/abnormal physiological events in each patient 
represented by said group; 

b) dividing each individual data into packets; 

c) transmitting the data packets of one or more patients to a remote site 
over a data network; and 

d) receiving said data packets in said remote site, processing the received 
data packets by a computerized system and displaying a representation of 
abnormal physiological events as a function of time and/or of changes in 
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normal physiological states and /or different stages of normal 
physiological states of each patient at said remote site. 

8. A method according to claim 7, further comprising: 

a) transmitting a request for one or more segments of the raw data streams 
that correspond to abnormal physiological events from the remote site to 
the data storage over the data network; and 

b) retrieving the requested segments from said data storage and displaying 
said requested segments at said remote site. 

9. A method according to claim 8, further comprising: 

a) transmitting instructions for a medical treatment to be given to the 
patient for terminating the abnormal physiological events from the 
remote site to a local computerized system located at the vicinity of 
said patient, over the data network; 

b) at said remote site, receiviag data packets that correspond to an 
individual group of data clusters that corresponds to physiological 
events after applying said medical treatment to said patient, 
processing the received data packets by the computerized system and 
displaying a representation of physiological events, obtained in 
response to said medical treatment, as a function of time and/or of 
changes in normal physiological states and /or different stages of 
normal physiological states of each patient; and 

c) optionally repeating the preceding steps, until abnormal physiological 
events are eliminated or substantially reduced. 

10. A method according to claim 1, wherein the raw data streams are received 
from one or more dry electrodes attached to the skin of the patient and having 
movement detectors for identifying time periods during which segments of raw 
data streams are associated with patient movements. 
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11. A method according to claim 1 
are selected from the group: 

- sleep apnea; 

- cardiac arrhythmia; 

- Ventricular Tachycardia; 

- Ventricular FibriUation/Flutter; 

- Epilepsy. 

12. A method according to claim 1, wherein the raw data streams are processed 
by performing the following steps: 

a) inputting said raw data streams to a cellular device having processing 
means being capable of identifying irregular physiological events of the 
patient; 

b) generating a characteristic trace that corresponds to regular physiological 
events of said patient; 

c) identifying selected raw data segments that comprise irregular 
physiological events; 

d) transmitting said selected raw data segments over a cellular network to 
another computational device; and 

e) continuing processing said selected raw data segments streams and 
representing the processed raw data by data clusters, each of which 
consisting of a set of unique properties that correspond to a different type 
of biomedical signal and to a different normal physiological state of said 
patient, by said another computational device. 

13. A method according to claim 12, further comprising: 

a) representing the abnormal physiological events and/or the individual 
attributes of each patient by an individual group of updated data clusters; 

b) dividing each individual group into data packets; 

c) transmitting the data packets of one or more patients to a remote site 
over a data network; and 
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, wherein the abnormal physiological events 
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d) receiving said data packets in said remote site, processing the received 
data packets by a computerized system and displaying a representation of 
abnormal physiological events as a function of time and/or of changes in 
normal physiological states and /or diEferent stages of normal 
physiological states of each patient at said remote site. 

14. A method according to claim 13, further comprising: 

a) for each patient, identifying, by another computational device, irregular 
physiological events; 

b) retrieving segments of raw data that corresponds to said irregular 
physiological events; 

c) generating a data set that comprises the retrieved segments and 
parameters indicating the irregularity of the irregular physiological 
events within said retrieved segments; 

d) transmitting said data set to a cellular device having processing means 
being capable of displaying said parameters and said retrieved segments, 
over a cellular data network; and 

e) displaying said parameters and/or said retrieved segments, on the 
display of said cellular device. 

15. A method according to claim 1, further comprising: 

a) Generating, in real-time, alert indications representing abnormal 
physiological events; and 

b) Automatically transmitting said indications to a physician over a data 
network, for determining/providing the required medical treatment. 

16. A method according to claim 1, further comprising dynamically controlling 
the operation of medical apparatus used for providing medical treatment to a 
patient being in an abnormal state in response to identified abnormal 
physiological events in said patient. 
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17. A method according to claim 16, wherein, the medical apparatus is selected 
from the group: 

- controllable drug dosage devices; 

- controllable gas/air delivery devices; 

- CPAP. 

18. A method according to claim 7, comprising transmitting instructions for 
dynamically controlling medical equipment, providing medical treatment to the 
patient, for terminating the abnormal physiological events from the remote site 
to a local computerized system located at the vicinity of said medical equipment, 

over the data network. 

19. A method according to claim 1, wherein the analyzed biomedical signal is an 
EGG signal, and the analysis of said signal comprises: 

a) acquiring the patient's raw EGG signal, by utilizing a corresponding 

electrode; 

b) performing adaptive segmentation of said raw EGG signal; 

c) extracting features from said raw EGG signal; 

d) clustering temporal and shape attributes of said features; and 

e) performing medical interpretation of the clusters. 

20. A method according to claim 19, wherein the adaptive segmentation 
comprises the steps of: 

a) generating several signals from the raw EGG signal, whereia in each 
signal a different frequency content being emphasized; 

b) summating the corresponding absolute values of said signals; 

c) filtering the resulting summation; and 

d) employing local maxima detection method, for identifying the R-peaks, P- 
peaks and T-peaks in said filtered resulting summation, said R, P and T 
peaks being utilized for characterizing the corresponding heartBeats 
Under Test (BUTs), said P and T peaks being utilized also for segmenting 
heartbeats. 
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21. A method according to claims 19 and 20, wherein the features extraction 
comprises the steps of: 

a) detecting 'R-R' time-intervals between each two consecutive R-peaks; 
and 

b) identifying 'P', 'Q', 'S' and T points of the corresponding BUTs, by 
utilizing said R-R' time-intervals, wherein said points characterize 
said BUT. 

22. A method according to claims 19 and 20, wherein the temporal shape, 
clustering and medicEd interpretation comprise the steps of: 

a) comparing the shape of the BUT to the shape of a Universal HeartBeat 
(UHB) shape; 

b) whenever the resemblance of said BUT to said UHB is smaller than a 
first Normalized Correlation threshold value NCjjji, ignoring said BUT; 
otherwise 

c) storing said BUT in a storage array; 

d) whenever required, averaging the stored BUTs, for obtaining a Patient- 
Specific Universal HeartBeat (PSUHB); 

e) replacing said UHB by said PSUHB; 

f) for each BUT, repeating steps a) to e), while utilizing a second 
Normahzed Correlation (NC) threshold value NCjy2' said second NC 
threshold value being larger than said first NC threshold value; 

g) whenever the resemblance of said BUT to the PSUHB is larger than the 
second Normalized Correlation threshold xalvieNC-pff2, storing said BUT as 
Normal heartbeat in a corresponding storage array; 

h) whenever required, averaguig the stored Normal heartbeats, for 
obtaining a Patient-Specific Normal HeartBeat (PSNHB); 

i) whenever the resemblance of said BUT is smaller than the second 
Normahzed Correlation threshold valne NCjjj2 > carrying out a Movement 
Decision Process (MDP); 
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j) if a corresponding movement is detected, by utilizing said MDP, 
discarding said BUT; otherwise 

k) comparing the shape of the BUT to the shape of a first aberrant 
Universal Pathological HeartBeat {UPBlgEp) shape and to a second aberrant 
Universal Pathological HeartBeat {UPBlj^); 

1) whenever the resemblance of said BUT to said UPBlg^F' suid/or to said 
UPB2j^p, is larger than a third Normahzed Correlation threshold 
value NCfHs, storing said BUT as Pathological heartbeat, the type of which is 
determined according to said UPBlgEp UPBlgj^p-, whichever yields a larger 
resemblance; otherwise 

m) if said BUT is found to be noise, discarding said BUT; otherwise 
n) whenever the resemblance of said BUT to said UPBlj^p, and/or to said 
UPBl^, is larger than a forth Normahzed Correlation threshold 
value #07^4, storing said BUT as Pathological heartbeat, the type of which is 
determined according to said UPBlj^p or UPBlj^p, whichever yields a larger 
resemblance; otherwise 

0) classifying said BUT as having unknown shape; 

p) whenever required, and for each type of Pathological he£irtbeat, 
averaging the stored Pathological HeartBeats, for obtaining a corresponding 
first Patient-Specific Universal Pathological HeartBeat (PSUPHBl) and a 
second Patient-Specific Universal Pathological HeartBeat (PSUPHB2); and 
q) for each additional BUT, repeating steps 1) to p), while replacing said 
UPBljuEF and UPBlj^Ephy said PSUPHBl and PSUPHB2. 

23. A method according to claim 22, in which the stored aberrant 
pathological BUTs are utiUzed for determining pathological sequences. 

24. A method according to claim 23, wherein the pathological sequence is 
selected firom the group: 

- VT; 
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- SVT; 

Bigeminy; 
Trigeminy 1; 
Trigeminy2; 
Tachycardia 
Bradycardia 

25. A method according to claim 22, wherein the BUTs are compared to the 
relevant shape by utilizing a corresponding Matched Filter, the parameters of 
which are updated whenever a new average Normal heartbeat or new 
Pathological heartbeat are calculated. 

26. A method according to claim 22, wherein identifying the R-peaks in the 
acquired ECG signal is carried out by utilizing the Wavelet Transform Process 
and several scales associated with said R-peaks. 

27. A method according to claim 1, in which the analyzed biomedical signal is 
an EEG signal, and the analysis of said signal comprises the steps of: 

a) acquiring the patient's raw EEG signal; 

b) Filtering the EEG signal, EGG and EMG signals; 

c) Adaptively segmenting the resulting filtered EEG and EOG signals into 
Quasi-Stationary Segments (QSSs); 

d) Extracting the required features from the resulting quasi-stationary EEG 
and EOG segments; 

e) Performing preliminary classification of each EEG QSS; 

f) Whenever required, detecting special temporal events contained within 
the EEG QSS; 

g) Performing special events classification, by using environment 
characteristics; 

h) If an EOG signal exists, detecting EOG REM, by employing related 
process; and 

i) Determining sleep staging. 
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28. A method according to claim 27, further comprising: 
a) performing a real-time analysis, by: 

a. 1) adaptive segmentation of said EEG signal, by utilizing a first 
and a second time-windows, the time-width and relative location being varied 
until a decision, regarding the optimized location of the boundaries of each EEG 
segment, is determined, based on the compaiison between the statistical 
properties of a jfirst EEG signal portion contained in the first time-window to 
the statistical properties of a second EEG signal portion contained in the second 
time-window; 

a.2) for each identified quasi-stationary segment, associating said 
segment to a High Frequency Group (HFG), or to a Mixed Frequency Group 
(MFG) or to a Low Frequency Group (LFG), by identifying its temporal and 
fi-equency characteristics and by employing fuzzy logic rules, said segments 
associated with same group form a cluster; 

a.3) whenever a segment is associated with HFG or LFG, analyzing 
its adjacent right-hand side and left-hand side segments, for identifying special 
events; 

a.4) whenever said segment significantly differs fi-om said adjacent 
segments, classifying said segment as a special event, which does not belong to 
the ongoing EEG signal; 

a.5) after first predetermined time-interval, averaging each one of 
the three related clusters; 

a.6) determming the sleep stage, being associated with said time- 
interval, by employing a set of R&K rules on the averaged clusters and on the 
detected special events; and 

a.7) repeating steps a.l) to a.6) for each segment. 

b) Performing an off-fine analysis, during which the system converges to the 
inherent sleep behavior of the monitored person, said analysis being employed 
for a second time-interval, the duration of which being larger than the duration 
of said first time-interv£d, comprising: 
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b.l) Performing initial division of the frequency groups into subgroups; 

b.2) employing the Baum-Welch process on the subgroups, for allowing 
learning the monitored person by averaging the corresponding subgroups; 

b.3) utilizing the averaged subgroups for refining the special event 
classification; 

b.4) detecting EMG energy lows using hierarchical fuzzy clustering 

process; 

b.5) utilizing the averaged subgroups for refining sleep staging process; 

and 

b.6) refining arousal detection by utilizing the results of the preceding 

steps. 

29. A method according to claim 1, in which the analyzed biomedical signals are 
Respiratory signals, and the analysis of said signal comprises the steps of: 

a) Detection and Smoothing of the Envelope (DSE) of the Chest effort 
signal, said effort signal, being referred to as the Master event; 

b) Peak Detection and Maximum Setting (PDMS) of the envelope signal; 

c) Identifying two consecutive global maxima points, said points defining 
the temporal boundaries of a corresponding segment of said envelope; 

d) Calculating the probability of a decrease of more than 20% in the 
magnitude of said segment; 

e) If a decrease of more than 20% is detected, for duration larger than 10 
seconds, employing a Fuzzy Logic Process, for evaluating the 
relevance of other events to said decrease, after which a decision is 
made regarding the probability of respiratory events; otherwise 

{) Identifying the trend of said segment, by assigning a first parabola 
function that mathematically approximates said segment and 
analyzing the characteristics of said parabola; 

g) If said trend is Downwards, assigning a second parabola function to 
the next detected segment; 
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h.) If the trend of said second parabola is Upwards, employing a process 
for adjoining the corresponding segments, said adjoined segments 
for min g a grand segment; and 

i) Employing step e) on the grand segment. 

30. A method according to claim 29, in which the evaluation process further 
comprises: 

a) associating a trapezoid-shaped Probability- Wise Window (PWW) to a 
corresponding event, the shape, time width and location of said PWW, 
in relation to the master event, being predetermined according to 
clinical criteria; 

b) utilizing said PWW for assigning a probability value to the 
corresponding event, according to the temporal location of the event in 
relation to the PWW associated with it; and 

c) calculating the probabihty of real respiratory events detection, based 
on the probabihty values being assigned to corresponding events. 

31. A method according to claim 29, in which the other events are the Abdominal 
Effort, Limb Movement and Oxygen Saturation. 

32. A method according to claim 29, in which the master event is selected from 
the group of: Abdominal Effort, Limb Movement, Oxygen Saturation and 
Chest Effort, to said remaining events are assigned a corresponding PWW, 
for allowing relating each one of said remaining events to said master event. 

33. System for momtoring normal/abnormal physiological events in one or more 
patients by the analysis of biomedical signals of each patient, comprising: 

a) One or more sensors, attached or allocated to said patient, being 
capable of generating and transmitting biomedical signals; 

b) a data storage for aggregating transmitted raw data streams, each of 
which representing a different type of biomedical signals of said 
patient; and 
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c) processing means, being in data communication with said sensors 
and/or said data storage, for performing adaptive segmentation of said 
raw data streams of said patient and assigning, to each segment, 
individual attributes being represented by attribute values 
determining an attribute domain, in which each segment being 
represented by a point, for generating a set of clusters in said 
attribute domain for each type of biomedical signal or for any 
combination thereof, such that each cluster consists of a combination 
of points determined by their relative location to other points, by 
assigning a set of property values to each point, and such that each 
property value corresponds to the degree of association of said point 
with one of the clusters, for obtaining the probabihty of the order of 
appearance of each point, in time, according to its property value and 
updating the property value of each point in each cluster using said 
probabihty, for updating each cluster by iteratively varying the 
combination of points included in each cluster according to their most 
updated property value and by including points derived from adaptive 
segmentation of further aggregated raw data streams of said patient, 
for associating at least one updated cluster with a normal/abnormal 
physiological state of said patient by using former knowledge, 
regarding normal/abnormal physiological states of said patient and/or 
of a reference group of patients, that is represented as reference 
clusters in said domain and for individually characterizing said 
patient by identifying normal/abnormal physiological states, 
associated with one or more updated clusters, and obtaining the 
probabihty of a change between normal/abnormal physiological states 
using said probabihty of the order of appearance. 



System according to claim 33, further comprising display means for 
displaying different stages of normal physiological states of the patient as a 
function of time. 
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35. System according to claim 33, in which the types of biomedical signals are 
selected from the group: 

- EEG signals; 

- EOG signals; 

- EGG signals; 

- Respiratory signals; 

36. System according to claim 33, further comprising: 

a) Means for dividing data dusters and/or the representation of their 
corresponding physiological states into data packets; 

b) means for transmitting said data packets of one or more patients to a 
remote site over a data network; 

c) means for receiving said data packets in said remote site; and 

d) means for processing the received data packets and displaying a 
representation of abnormal physiological events as a function of time 
and/or of changes in normal physiological states and/or different 
stages of normal physiological states of each patient at said remote 
site. 

37. System according to claim 36, further comprising: 

a) means for transmitting a request for one or more segments of the raw 
data streams that correspond to abnormal physiological events from 
the remote site to the data storage over the data network; and 

b) means for retrieving the requested segments from said data storage 
and displaying said requested segments at said remote site. 

38. System according to claim 33, further comprising one or more dry electrodes 
attached to the skin of the patient and having movement detectors for 

generating and forwarding raw data streams. 



39. System according to claim 33, further comprising: 
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a) a computational device having processing means for identifying 
abnormal physiological events in a patient, and being capable of 
generating a data set that comprises retrieved segments and 
parameters indicating the irregularity of the irregular physiological 
events within said retrieved segments and of transmitting said data 
set; and 

b) a cellular device having processing means being capable of receiving 
said data set over a cellular data network, and a display for displaying 
said parameters and said retrieved segments. 

40. A method according to claim 1, wherein the former knowledge is capable of 
being continuously updated by manually modifying the characterization 
results. 

41. A method according to claim 28, wherein the relative energy of the high 
frequency group, mixed frequency group and low frequency group, is 
displayed as a fiinction of time, in the form of two-dimensional graph, the 
graph of each one of the frequency groups being displayed by different color, 
for allowing an observer distinguishing between said frequency groups. 

42. System for allowing a plurality of distributed users retrieving and 
displaying, essentially simultaneously, data related to preferred aspects of 
physiological condition of a plurality of distributed patient's data records or 
controlling the monitor/treatment device, comprising: 

a) at least two one-Lab-Server systems, each one-Lab-Server comprising: 
a.l) at least one gateway, to which at least one monitor and/or at least 
one treatment device are connected, each one of said at least one monitor 
and at least one treatment device may optionally be connected to at least 
one patient, via corresponding electrodes, or sensors, said gateway being 
capable of preprocessing raw data representing the physiological signals 
of said patient and beiag monitored by said monitor, and formatting the 
preprocessed raw data, so as to allow said data to be forwarded to a 
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remote user who requests said formatted data, said gateway being 
capable of forwarding control signals to said treatment device, for 
improving the monitored physiological condition of said patient; 
a.2) at least one user workstation, for allowing a physician access to the 
medical information contained within said system, after checking the 
authority level of said physician, and forwarding control signals, via an 
intermediator, and via corresponding gateway, to corresponding 
treatment apparatus connected thereto; and 

a.3) A Lab Server, being said intermediator, and being capable of 
estabHshing bi-directional communication with said at least one gateway 
and with said at least one user workstation, via a data network, said 
formatted data is forwarded to said Lab Server, from a corresponding 
Gateway and via said network, said Lab Server being programmed to 
process and analyze said formatted data, said at least one user 
workstation being capable of forwarding to said Lab Server, via said 
network, at least one request and control signals, said request being 
utilized for retrieving, from, or via, said Lab Server, and displaying, on 
said user workstation, processed and analyzed data related to specific 
aspects of the physiological condition of a specific patient, said specific 
aspects being defined by a user of said at least one user workstation, and 
said control signals being relayed, by said Lab Server, to corresponding 
treatment device, for affecting the treatment process applied to a patient 
connected thereto, said Lab Server comprising a database for storing said 
raw, formatted and processed/analyzed data; and 
b) A central Server, being capable of establishing bi-directional 
communication with the Lab Server of each one of said one-Lab-Server 
systems, via the data network, and with each one of the user work 
stations, said user work stations being capable of forwarding to said 
central server login information and request to receive specific medical, or 
physiological, data of specific patient, said central server being 
programmed to identify, using a central database contained therein, the 
corresponding Lab Server which contains said specific medical data of 
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specific patient, after which said central server routes, or relays, said 
workstation to said Lab Server for allowing said workstation to obtain 
said specific medical data of specific patient. 

. System according to claim 46, wherein the One-lab server system is self- 
contained. 
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